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Chapter  1 
Introduction 


The  problem  of  designing  a  feedback  system  to  satisfy  a  performance  specifica¬ 
tion  and  to  be  robustly  stable  against  plant  modelling  errors  is  of  fundamental 
importance  in  control  theory.  Indeed,  over  the  last  10  yean,  this  problem  has 
received  renewed  attention,  particularly  in  the  context  of  multivariable  systems 
(e.g.  [D0S8I],  [SLH81],  [Fra82],  [IYa87],  [Doy87],  [Fra87]).  Singular  value  analy> 
sis,  baaed  upon  the  singular  value  norm  of  a  matrix  transfer  function,  has  served 
aa  an  important  tool  for  characterising  both  performance  and  stability  robustness 
properties  of  a  multivariable  feedback  system,  and  has  allowed  the  generalisa¬ 
tion  of  many  useful  concepts  from  classical  feedback  theory  for  scalar  systems  to 
multivariable  systems. 

Unfortunately,  singular  value  analysis  techniques  are  not  applicable  to  many 
design  problems  of  interest.  For  example,  singular  value  robustness  analysis  is 
moot  useful  when  uncertainty  is  isolated  at  only  one  point  in  the  system.  When 
uncertainty  and  modelling  errors  are  present  at  several  points  in  the  system, 
singular  value  analysis  tends  to  yield  estimates  of  stability  robustness  that  are 
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either  optimistic  (by  testing  robustness  against  only  one  modelling  error  at  a 
time)  or  pessimistic  (by  testing  against  a  broader  class  of  uncertainty  than  is 
actually  present).  Furthermore,  singular  value  analysis  is  useful  only  for  test¬ 
ing  nominal  performance  properties  of  a  system;  obviously,  performance  should 
also  be  robust  against  modelling  error.  It  may  be  shown  [DWS82]  that  the  ro¬ 
bust  performance  problem  is  equivalent  to  one  of  stability  robustness  with  an 
additional  source  of  uncertainty  that  represents  the  performance  specification. 
Hence,  robust  performance  problems  reduce  to  problems  of  stability  robustness 
with  respect  to  several  sources  of  uncertainty,  and  the  deficiencies  of  singular 
value  analysis  described  above  also  apply  to  the  analysis  of  robust  performance. 

To  address  the  limitations  of  singular  value  analysis,  Doyle  ([Dqy82],  [DWS82]) 
introduced  the  structured  singular  value,  an  analysis  tool  that  directly  addresses 
the  problem  of  stability  robustness  against  several  sources  of  modelling  uncer¬ 
tainty,  and  thus  also  addresses  the  robust  performance  problem.  Essentially, 
the  structured  singular  value  provides  a  precise  stability  margin  against  several 
simultaneous  sources  of  modelling  error. 

With  the  aid  of  singular  value  and  structured  singular  value  analysis,  several 
methodologies  for  the  design  of  robust  multivariable  feedback  systems  have  been 
proposed.  Among  these  are  LQG/LTR  ([D0S8I],  [StA87]),  H°°  optimisation 
[Fra87],  and  structured  singular  value  synthesis  [Doy85].  At  present,  each  of  these 
methodologies  suffers  from  at  least  one  shortcoming.  For  example,  the  LQG/LTR 
methodology  is  applicable  only  to  robustness  problems  with  one  source  of  mod¬ 
elling  error.  This  shortcoming  is  overcome  by  the  structured  singular  value  syn¬ 
thesis  methodology;  however,  the  latter  methodology  is  experimental  and  tends 
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to  generate  exceedingly  high  order  compensators. 

It  seems  clear  that  additional  insights  into  the  properties  of  multivariable 
systems  with  several  sources  of  uncertainty  will  be  needed  before  a  completely 
satisfactory  design  methodology  can  be  developed.  In  particular,  multivariable 
systems  can  possess  robustness  difficulties  having  no  analogue  in  their  scalar 
counterparts.  For  example,  it  has  been  observed  by  several  researchers  ([Ste84], 
[Ste85],  [Ste87],  [Doy87],  [NeM87],  [SMD88],  [Fre89a],  Fre89b})  that  systems  with 
ill-conditioned  plants  can  cause  robustness  difficulties  when  modelling  error  is 
present  at  more  than  one  loop  location.  Ill-conditioned  plants  arise  in  fields  as 
diverse  as  aircraft  control  [Enn87]  and  chemical  process  control  [SMD88].  From 
a  systems  viewpoint,  ill-conditioning  at  a  frequency  means  that  the  gun  of  the 
plant  exhibits  a  strong  directional  dependence;  i.e.,  certain  input  sinusoids  at 
a  frequency  for  which  the  plant  is  ill-conditioned  will  be  amplified  much  more 
than  will  others.  A  thorough  understanding  of  robustness  problems  caused  by 
plant  ill-conditioning  is  not  yet  available,  nor  is  a  design  methodology  capable  of 
coping  with  these  difficulties. 

The  work  performed  under  this  contract  represents  substantial  progress  to¬ 
ward  understanding  the  robustness  properties  of  multivariable  systems  with  sev¬ 
eral  sources  of  uncertainty,  and  toward  design  problems  associated  with  an  ill- 
conditioned  plant.  The  springboard  for  the  present  work  is  provided  by  the 
results  of  [Fre89a,  and  [Fre89b].  To  summarise  these  results  briefly,  the  stabil¬ 
ity  robustness  problem  with  two  sources  of  uncertainty  is  analysed  to  determine 
when  a  system  that  is  robust  against  each  source  of  modelling  error  (assuming 
that  the  other  is  not  present)  can  nevertheless  be  destabilised  by  simul- 


3 


taneous  modelling  errors.  These  results  are  then  applied  to  the  ill-conditioned 
plant  problem.  A  sensitivity  analysis  is  performed  to  detect  small  modelling 
errors  that  strongly  affect  system  stability  robustness  and  performance.  These 
modelling  errors  are  then  used  to  derive  bounds  upon  the  stability  robustness 
margin,  as  quantified  by  the  structured  singular  value.  The  significance  of  these 
bounds  is  that  they  express  the  interrelation  among  the  plant  directionality  prop¬ 
erties,  the  design  specifications,  and  the  compensator.  This  information  suggests 
a  strategy  for  compensator  design  to  achieve  robust  performance  and  stability 
despite  plant  ill-conditioning. 

The  research  described  in  this  report  contributes  to  the  understanding  of 
multivariable  robustness  problems  in  several  ways.  Our  overall  approach  is  sum¬ 
marised  as  follows.  We  shall  study  robustness  problems  that  may  be  analysed 
using  the  Doyle’s  structured  singular  value  (SSV).  Because  calculation  of  the 
SSV  involves  numerical  optimisation  procedures,  it  is  difficult  to  obtain  insight 
into  the  relation  between  its  sise  and  the  plant,  compensator,  and  design  spec¬ 
ifications.  Such  insight  is  useful  for  several  reasons.  First,  one  would  like  to 
identify  combinations  of  plant  properties  and  design  specifications  that  are  in¬ 
herently  difficult  to  satisfy.  Second,  one  would  like  to  have  the  capability  to 
manipulate  feedback  properties,  such  as  robustness  and  performance,  by  modify¬ 
ing  the  compensator.  Our  approach  to  obtaining  such  insight  is  to  derive  bounds 
and  approximations  to  the  SSV.  To  be  useful,  these  bounds  and  approximations 
must  be  both  reasonably  accurate,  and  should  also  display  the  dependence  of 
the  SSV  upon  plant,  compensator,  and  design  specifications.  Once  such  bounds 
and  approximations  have  been  derived,  they  are  then  analysed  to  provide  insight 
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into  fundamental  design  limitations,  and  to  suggest  compensation  strategies.  In 
particular,  we  incorporate  the  information  gleaned  from  the  bounds  into  a  se¬ 
lection  procedure  for  the  weightings  used  in  LQG/LTR,  H 00 ,  and  H*  synthesis 
procedures. 

The  remainder  of  this  report  is  outlined  as  follows.  In  Chapter  2  we  shall 
derive  bounds  upon  the  SSV  for  robustness  problems  with  two  or  three  sources 
of  uncertainty.  The  bounds  with  two  sources  of  uncertainty  improve  those  ap¬ 
pearing  in  [Fre89a],  and  [Pre89b],  while  those  with  three  sources  of  uncertainty 
are  new.  In  Chapter  3,  we  analyse  a  specific  robustness  problem,  that  of  output 
performance  with  respect  to  input  uncertainty  which  is  diagonal  or  block  diago¬ 
nal.  Using  a  sensitivity  analysis,  as  well  as  the  bounds  from  Chapter  2  for  three 
sources  of  uncertainty,  we  study  the  design  difficulties  posed  by  an  ill-conditioned 
plant  for  such  a  robustness  problem.  In  Chapter  4  we  develop  a  framework  for 
analysing  robustness  with  respect  to  several  sources  of  uncertainty.  We  identify 
certain  interaction  parameters, ,  whose  sise  determines  the  extent  to  which  inter¬ 
actions  among  several  modelling  errors  can  cause  robustness  difficulty,  and  use 
these  parameters  to  derive  bounds  upon  the  SSV.  Chapters  5  and  6  are  devoted 
to  methods  of  incorporating  insights  obtained  from  our  analyses  into  the  JET00 
and  LQG/LTR  design  methodologies.  Chapter  7  contains  directions  for  further 
research. 
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Chapter  2 

Bounds  on  the  Structured 
Singular  Value  with  Two  or 
Three  Sources  of  Uncertainty 


2.1  Introduction 

The  purpose  of  this  chapter  is  to  present  bounds  upon  the  structured  singular 
value  (SSV)  with  two  or  three  sources  of  uncertainty.  For  two  sources  of  uncer¬ 
tainty,  our  bounds  improve  those  in  [Fre89a]  and  [Fre80b].  For  three  sources  of 
uncertainty,  our  bounds  can  be  shown  to  improve  those  of  Deznmel  [Dem88].  Ap¬ 
plications  of  our  bounds  to  specific  classes  of  robustness  problems  will  be  found 
in  Chapter  3  of  this  report. 

The  remainder  of  this  chapter  is  organised  as  follows.  Section  2.2  briefly 
reviews  the  structured  singular  value  and  its  properties.  Section  2.3  contains  our 
main  results,  including  both  lower  and  upper  bounds  upon  the  3-block  structured 
singular  value.  These  bounds  are  expressed  in  terms  of  a  set  of  interaction 
parameters  that  essentially  determine  how  several  uncertainties  can  interact  to 
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cause  robustness  difficulties.  The  research  described  in  this  chapter  is  discussed 
in  [ChFSOc]  . 

2.2  Properties  of  the  Structured  Singular  Value 

We  shall  now  briefly  review  those  properties  of  the  structured  wingwUr  value  used 
in  this  chapter;  a  complete  discussion  and  examples  are  found  in  [Doy82].  The 
first  step  in  the  structured  singular  value  analysis  is  to  rearrange  the  block  dia¬ 
gram  the  feedback  system  into  the  form  shown  in  Figure  2.1.  The  uncertainty 
matrix  A(«)  is  assumed  to  lie  in  the  set 


A  =  {  A(a)  :  A(«)  =  diag[A1(«),  Aj(s),  •  •  • ,  A*(«)],  A <  €  C4****  and  stable  } 

(2-2-1) 

It  is  often  convenient  to  assume  that  system  uncertainties  have  been  scaled  to 
satisfy  the  upper  bound  <  7,  Vw,  Vt,  and  to  introduce  the  set  A7  = 

{  A(s) :  A(s)  €  A  and  <T[A(jw)j  <  7,  Vw  }.  The  interconnection  matrix 

Af(.)  =  [**/(•)]  =  «,(*)  e  C*’"!!  (2J.2) 

is  stable  if  the  feedback  system  is  nominally  stable.  At  each  frequency,  the 
structured  singular  value,  denoted  p[A f],  is  defined1  by  [Doy82],  [DWS82])* 

u[M]  _  f  0  det[I  +  MAj  *  OVA  €  A 

^  \  l/(minA6A{^[A] :  det\I  +  M A]  =  0})  otherwise 

(2.2.3) 

1  Throughout  this  sots,  «i  suppress  dependence  upon  frequency  whenever  ^pnpriiU. 
aBy  a  mild  abuse  of  notation,  tbe  symbol  A  will  occasionally  be  used  to  denote  the  set  of 
constant  complex  matrices  with  block  diagonal  structure  (2.2.1). 
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Hence,  by  definition,  the  reciprocal  of  the  structured  singular  value  is  a  fre¬ 
quency  dependent  stability  margin  with  respect  to  uncertainty  of  the  form  (2.2.1), 
and  a  nominally  stable  feedback  system  will  be  robustly  stable  against  all  uncer¬ 
tainty  in  the  set  A,  if  and  only  if  n[M(ju)]  <  I/7,  Vw. 

The  value  of  n[M]  can  be  calculated  as  follows.  Define  the  set 

P  =  {  diag[dilklt  da/*,,  •  •  • ,  d*J*J  :  4  €  R,  dt  >  0  }  (2.2.4) 

where  /*,  denotes  the  identity  matrix  of  dimension  He*.  It  was  shown  in  [Doy82] 
that 

H[M]  <  wf>  9[DMD~l]  (2.2.5) 

In  particular,  the  upper  bound  infnepff[DMD_1]  is  equal  to  p[M]  when  the 
system  has  three  or  fewer  blocks  of  uncertainty,  namely,  k  <  3.  This  fact  may 
be  used  to  derive  bounds  upon  the  structured  singular  value. 


2.3  Lower  and  Upper  Bounds  of  /i[Af] 


In  this  section  we  present  both  lower  and  upper  bounds  on  the  3-block  structured 

singular  value.  For  convenience,  we  shall  first  denote  ay  —  and  introduce 

the  following  parameters  that  will  be  used  subsequently  in  deriving  the  bounds. 

m(»\j)  =  y/oijofii  Vi/; 

J»(«. /,*)  =  {oijOjkOu)111  Vi  /;'  /  fc 

Hi  =  max*  an  (2.3.1) 

M*  =ma Xijjtfaj) 

Hi  =  mAXiJikH[i,j,k) 

Interpretations  of  these  parameters,  which  we  shall  term  ‘interaction  parame¬ 
ters,"  will  become  available  as  we  develop  bounds  upon  the  structured  singular 
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value.  For  motivation,  let  us  postulate  an  analysis  procedure  wherein  we  test 
robustness  against  uncertainties  taken  first  one  at  a  time,  then  two  at  a  time 
and  finally  all  three  at  a  time.  We  wish,  in  particular,  to  understand  when  two 
(or  three)  uncertainties  can  interact  to  cause  robustness  difficulties  even  though 
robustness  is  good  with  respect  to  uncertainties  taken  one  (or  two)  at  a  time. 
It  is  well-known  that  the  stability  margin  against  A <  alone  is  equal  to  l/o«  ( 
[DWS82].  Hence  the  parameter  Hi  measures  stability  robustness  against  the  un¬ 
certainties  taken  one  at  a  time.  Suppose  next  that  we  consider  the  effects  of 
uncertainties  taken  two  at  a  time.  With  no  loss  of  generality,  we  may  study  this 
problem  by  setting  k  —  2  in  (2.2.1-2). 

Proposition  2.8.1  Let  k  =  2  m  (8.8.1-8),  and  suppose  that  the  infimum  in 
(8.8.5)  is  achieved.  Then 

max  {  nu  Mi }  <  v\M\  <  Mi  +  Mi  (2.3.2) 

Proof:  In  order  to  prove  the  lower  bound,  we  note  that  #*[Af]  =  for 

some  t>  -  diagldilk^dilkt]  €  P.  It  then  follows  that  p[Af]  >  (di/dj)oiJt  Vi,j. 
Hence  h[M]  >  <rM,  and  p*[M]  >  (dj/d,-)oy  '{dj/d^o#  =  (*>/).  thus  establishing 
the  lower  bound.  The  proof  for  the  upper  bound  is  found  in  [FYe89a]. 


Remark  2.8.1:  The  upper  bound  in  (2.3.2)  first  appeared  in  [FreSOa],  and 
the  lower  bound  is  a  tighter  version  of  that  in  [£Ve89a].  The  statement  and 
derivation  of  the  improved  lower  bound  are  due  to  C.  N.  Nett  (personal  commu¬ 
nication).  The  proof  of  this  improved  bound  is  significantly  simpler  than  that  of 
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the  bound  in  [Fre89a],  although  the  proof  in  [Fre89a]  is  conceptually  appealing 
in  that  it  provides  the  set  of  smallest  destabilising  uncertainties. 

Obviously,  Proposition  2.3.1  can  be  extended  to  analyse  the  effects  of  any 
two  uncertainties  in  the  general  Jk-block  structured  singular  value  problem.  The 
parameter  p(i,j)  essentially  determines  how  the  uncertainties  A<  and  Ay  can  in¬ 
teract  to  produce  robustness  difficulties;  hence  motivating  the  terminology  inter¬ 
action  parameter.  Note  in  particular  that  a  stable  system  whose  interconnection 

cannot  be  destabilised  by  either 
uncertainty  acting  alone;  as  pointed  out  in  [Fre89a]  and  [NeU88[,  p[M]  =  p%  in 
this  case. 

We  now  extend  this  result  to  the  3-block  structured  singular  value  problem. 
Using  Nett’s  technique  to  extend  the  lower  bound  in  (2.3.2)  is  straightforward. 

Proposition  2.8.2  Let  k  —  3  in  (2.2.1 ,2),  and  evppoee  that  the  infimum  in 
(2.2.5)  is  achieved.  Then, 

p[M]  >  max  {pu  pa  >  (2.3.3) 

■ 

Proof:  It  suffices  to  show  p[M]  >  ps.  Following  the  same  reasoning  as  in  the 
proof  of  Proposition  3.1  leads  to  p3[M]  >  {di/dj)oy  •  ( dj/dk)ojk  •  (d*/d*)o*,  for 
some  di,  dt,  and  d*.  Thus,  p*[M]  >  p*(i,j,k),  Vi  /  j  k.  m 

Remark  2.8.2:  Note  that  this  proof  technique  consists  of  finding  a  lower  bound 
upon  infflep  0\DMD~l],  and  may  be  extended  to  derive  similar  lower  bounds 
upon  this  quantity  for  the  case  of  k  >  3.  Of  course,  this  technique  results  in  a 
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matrix  has  the  form  M(s)= 


lower  bound  upon  p[M\  only  if  the  upper  bound  in  (2.2.5)  is  an  equality,  and  this 
is  guaranteed  only  for  k  =  2,3.  Most  algorithms  (e.  g.,  see  [FaT86],  [FaT88]) 
for  computing  the  structured  singular  value  actually  compute  this  upper  bound, 
which  Doyle  has  conjectured  is  within  15%  of  the  true  value  of  p\M]  ([Doy82], 
[DLP86]).  Hence,  if  Doyle’s  conjecture  is  correct,  one  might  argue  that  extensions 
of  Proposition  2.3.2  to  &  >  3  remain  useful. 

Remark  2.3.3:  Alternately,  one  can  also  apply  the  techniques  of  [Fre89a]  for 
the  case  k  =  2  to  derive  lower  bounds  for  k  >  3.  As  in  the  case  k  =  2,  the  lower 
bound  obtained  in  this  way  is  not  as  sharp  as  (2.3.2).  However,  the  technique  is 
in  principle  generalisable  to  obtain  lower  bounds  for  arbitrary  k. 

Next,  we  present  upper  bounds  upon  p[M).  Toward  this  end,  it  is  necessary 
to  consider  two  cases. 

Proposition  2.8.8  Letk  —  Z  in  (t.8.1,8). 

(a)  Suppose  Ms  >  0.  Then 

l*[M]  <  Mi  +  Ms  + 

Ms 

If  Mi/Ms  <  1,  then 

mM  <  Mi  +  Ms  +  Ms 

(b)  Suppose  ms  >  0.  Then 

M{Af]  <  Mi  +  Ms  +  V3  max  {ms,  (— )l,t  > 

Ms 

If  Mi/Ms  >  1.  then 

M(Af  ]  <  Mi  +  (1  +  V3  )ms 


(2.3.4) 

(2.3.5) 

(2.3.6) 

(2.3.7) 
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Proof  of  (a):  With  no  loss  of  generality,  we  wwnime  Mi  =  m(1»2,3).  Our  proof  pro¬ 


ceeds  with  the  decomposition  M  =  Mo+Ml+Ms,  where  Mo  =  diag[M\u  Mm,  Mm], 


'  0  Mis  O' 

o 

o 

Mx  = 

0  0  Mss 

and  Ms  = 

Mil  o  o  | 

Mn  0  0 

.  0  Ms,  0  J 

(2.3.8) 


Note  first  that 


Ms  =  rnf  o[DM\D~l\ 

.  t  /  d i  d,  d$  1 

“  ma  1  Zca'  Zan'  ) 

The  solution  to  this  minimisation  problem  is  obtained  by  setting  Ji  =  1,  d,  = 
ow/m(  1,2,3),  <4  =  m(1.2,3)/<7si,  leading  to  Ms  =  m(1.2,3).  Next,  notice  that 

m[M]  <  ini  9[DMD-1]  <  Mi  +  *[£Mi£-1]  +  *[£m,£_11 


where 


o[DMif)~l\  =  max  |  ^r<7u,  | 

[  M*(l,2)  m»(2,3)  m»(M)  1 
1  Ms  ’  Ms  *  Ms  J 


=  max 


<  A 


Ms 


This  proves  (2.3.4),  and  (2.3.5)  immediately  follows. 

Proof  of  (b):  Without  loss  of  generality,  we  wswnme  that  Ms  —  m(1*2).  Our 
proof  proceeds  with  the  decomposition  M  =  Mo  +  M,  +  M«,  where  Mo  = 


d»ag[Mn,  Mss,  Mss], 


0  Mi,  O' 

0  0  Mis 

M,= 

Mn  o  o 

and  M4  = 

0  0  Mm 

0  0  0 

Msi  Ms,  0 
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Clearly,  n[M)  <  Mi  +  infD&>  {o[DMtD-l\  +  9[DMAD~l]  >.  Setting  4  =  1, 4  = 
tfu/Mi  yield® 

M[-W]  <  Mi  +  Ms  +  W  ,  fr[DAf4.D-1] 

4iA“«ji4i»4 

To  solve  the  minimisation  problem  in  the  right  hand  side,  note  that 


inf  9[DM4D'1]  = 


^  max  |dgff[Msi,  (1/4)-Mm|*  ijMn  } 

<  inf  max  |  <k\Joh  +  (l/A)1^  *  (l/dajV^u +  «§*»} 


The  last  minimisation  problem  has  the  solution  A  =  [(?u  +  ^°u)/(ah  + 
(l/A)*oJ*)]1/4.  Substituting  the  value  of  A  yields 

al6M.fi-1]  <  ^*(2,3)  +  ».4(1,3)  +  +  <‘*(1^2)  j  ' 


<  w 


»,+(a3’' 


M?  +  2- 
Ms 


This  proves  (2.3.6),  and  (2.3.7)  follows  immediately. 


Just  as  the  parameter  Ms  measures  how  the  effects  of  two  uncertainties  can 
interact  to  cause  robustness  difficulties,  so  does  Ms  measure  the  problems  caused 
by  interactions  among  all  three  uncertainties.  Indeed,  suppose  that  the  system 
interconnection  matrix  has  the  form  of  Afi(s),  or  Afs(s),  defined  as  in  (2.3.8). 
Since  Mi  =  Ms  =  0,  neither  of  these  systems  can  be  destabilised  by  uncertainties 
acting  individually  or  in  pairs.  However,  each  system  can  be  destabilised  by 
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a  combination  of  Ax,  As,  and  As.  The  following  corollary  is  an  immediate 
consequence  of  (2.3.3)  and  (2.3.4). 
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Corollary  2.S.4  Let  Mi  and  Ms  be  defined  ae  in  (t.9.8).  Then , 

H 

m[Mi]  =  m(1,2,3) 


0>) 


/i[Ms]  =  p(l,3,2) 


(2.3.0) 


(2.3.10) 


Hence  the  first  upper  bound  (2.3.4)  is  useful  in  the  case  that  the  most  significant 
difficulty  is  due  to  the  interaction  of  all  three  uncertainties,  while  the  second 
upper  bound  (2.3.6)  is  useful  in  the  case  that  the  most  significant  difficulty  is 
due  to  the  interaction  of  just  two  of  the  uncertainties.  An  interesting  open 
problem  is  whether  a  single  upper  bound  exists  which  is  useful  in  both  cases  and 
is  not  prohibitively  messy. 


2.4  Concluding  Remarks 

In  this  chapter  we  have  derived  both  lower  and  upper  bounds  upon  the  struc¬ 
tured  singular  value  with  respect  to  three  blocks  of  uncertainty.  Our  bounds  are 
expressed  in  terms  of  a  set  of  parameters  that  determine  how  two  or  three  un¬ 
certainties  can  interact  to  cause  robustness  difficulties.  These  bounds  appear  to 
give  reasonably  tight  estimates  of  the  structured  singular  value.  Specifically,  the 
upper  bounds  were  shown  to  be  within  a  factor  of  three  of  the  lower  bound  on 
one  occasion,  while  within  a  factor  of  2  +  y/%  (»  3.73)  on  the  other.  It  is  interest- 
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ing  to  note  that  our  bounds  are  substantially  tighter,  and  have  more  interesting 
interpretations,  than  those  developed  in  [Dem88]. 

More  importantly,  our  bounds  may  be  used  to  study  various  robustness  analy¬ 
sis  and  design  problems.  In  applying  these  bounds,  we  found  that  the  blocks  My 
of  the  interconnection  matrix  will  be  mutually  interrelated,  and,  as  in  [Fre88a], 
it  is  important  to  analyse  this  interrelation  to  obtain  design  insight. 

We  conclude  this  chapter  by  pointing  out  the  potential  extensions  of  these 
bounds  to  the  lb-block  structured  singular  value,  with  k  >  3.  Many  of  the  present 
results  can  be  extended;  however,  the  technique  used  to  derive  the  lower  bounds 
may  not.  The  major  obstacle  in  achieving  this,  of  course,  is  the  invalidity  of  the 
equality  n[M\  —  inf  Dev  9[DMD~l ]  in  the  general  case. 
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Chapter  3 

Robust  Performance  with 
Respect  to  Diagonal  Input 
Uncertainty 


3.1  Introduction 

In  this  chapter  we  shall  apply  the  results  of  Chapter  2  to  a  specific  class  of 
robust  performance  problems.  For  motivation,  consider  the  problem  of  main¬ 
taining  the  output  sensitivity  function  small,  to  achieve  disturbance  attenuation, 
despite  the  presence  of  unstructured  multiplicative  uncertainty  at  the  plant  in¬ 
put.  Analysis  of  this  problem  has  exposed  potential  design  difficulties  when  the 
plant  transfer  function  matrix  is  ill-conditioned  at  frequencies  near  crossover.  By 
“ill-conditioned,*  we  mean  that  the  gain  of  the  plant,  at  a  frequency  of  interest, 
exhibits  a  strong  directional  dependence:  same  input  signals  will  be  amplified 
much  more  than  will  others.  Although  no  conclusive  proof  has  yet  appeared, 
anecdotal  evidence  suggests  that  ill-conditioned  plants  may  be  inherently  dif¬ 
ficult  to  robustly  control,  at  least  for  certain  types  of  plant  uncertainty,  and 
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that  the  sise  of  the  plant  condition  number  is  an  indicator  of  the  degree  of  the 
difficulty  in  achieving  robustness  [Fre89a],  [Fre89b],  [NeM87],  [SMD88],  [Ste85], 
[Ste87],  [SkM87], 


The  tentative  conclusion  that  an  ill-conditioned  plant  is  inherently  difficult 
to  robustly  control  depends  critically  upon  the  assumption  that  the  input  uncer¬ 
tainty  is  unstructured,  and  can  therefore  introduce  coupling  among  different  plant 
inputs.  Often,  however,  physical  considerations  dictate  that  it  is  more  reasonable 
to  assume  that  uncertainty  cannot  introduce  such  coupling,  or  can  introduce  cou¬ 
pling  only  among  a  subset  of  the  inputs.  Examples  arise  in  diverse  applications, 
such  as  process  control  [SkM87],  [SMD88]  and  aircraft  control  [Enn87].  This 
property  is  modelled  mathematically  by  assuming  that  the  input  uncertainty  has 
a  diagonal  or  block  diagonal  structure,  hi  (SMD8S],  it  is  demonstrated  via  ex¬ 
amples  that  ill-conditioned  plants  can  sometimes  cause  design  difficulties  when 
input  uncertainty  is  constrained  to  be  diagonal  (i.e.,  to  introduce  no  coupling 
between  inputs),  but  that  sometimes  no  such  difficulty  is  encountered.  Appeal¬ 
ing  physical  explanations  of  this  phenomenon  are  presented  in  [SMD88],  and  the 
extent  of  the  potential  robustness  difficulty  depends  upon  the  input  directional¬ 
ity  properties  of  the  plant,  as  well  as  its  condition  number.  It  is  also  remarked 
[SkM87]  that  the  plant  relative  gain  array  may  be  a  useful  indicator  of  potential 
robustness  difficulty  for  this  problem.  The  analysis  in  [SkM87],  [SMD88]  focuses 
in  particular  upon  robustness  difficulties  that  can  be  encountered  through  use  of 
a  compensator  that  explicitly  inverts  the  plant  model. 

Our  purpose  in  this  chapter  is  to  provide  a  framework  useful  for  analysing  the 
degree  of  difficulty  inherent  in  the  robust  performance  problem  when  the  input 
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uncertainty  is  modelled  as  a  diagonal  or  block  diagonal  matrix.  An  outline  of 
our  results  follows.  In  Section  3.2,  we  introduce  notation  and  define  the  robust 
performance  problem  of  maintaining  the  output  sensitivity  function  small  despite 
the  presence  of  block  diagonal  uncertainty  at  the  plant  input.  Section  3.3  con- 
tuns  initial  results:  we  calculate  the  first-  and  second-order  effects  of  the  input 
uncertainty  upon  the  output  sensitivity  function,  and  point  out  conditions  on  the 
closed-loop  shapes  which  must  be  satisfied  if  these  effects  are  to  be  kept  small.  In 
Section  3.4,  we  review  the  use  of  the  structured  singular  value  [Doy82]  as  a  test 
for  robust  performance.  Since  the  interrelation  among  the  plant,  compensator, 
uncertainty  structure,  and  sise  of  the  structured  singular  value  is  rather  compli¬ 
cated,  we  use  the  frequency  dependent  bounds  upon  the  structured  wingnlar  value 
derived  in  Chapter  2  to  provide  useful  insights  into  this  interrelation.  These 
bounds  are  stated  in  terms  of  interaction  parameters ;  keeping  these  parameters 
small  is  approximately  equivalent  to  satisfying  the  goal  of  robust  performance. 
Furthermore,  the  interaction  parameters  are  very  closely  related  to  the  first- 
and  second-order  effects  of  input  uncertainty  that  were  calculated  in  Section  3.3. 
Next,  in  Section  3.5,  we  analyse  the  problem  of  keeping  the  interaction  parame¬ 
ters  small,  and  show  that  the  difficulty  of  this  problem  depends  upon  a  robustness 
indicator  which  is  a  function  of  both  the  plant  condition  number  and  its  input 
directionality  properties.  In  Section  3.6,  we  briefly  discuss  and  compare  the  rel¬ 
ative  gun  array  to  our  new  robustness  indicator.  Examples  are  given  in  Section 
3.7,  and  concluding  remarks  are  made  in  Section  3.8.  The  research  described  in 
this  chapter  is  discussed  in  [ChF89a]. 
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3.2  The  Robust  Performance  Problem 


Consider  the  linear  time-invariant  feedback  system  depicted  in  Figure  3.1.  The 
transfer  functions  P(s)  and  F(s)  are  those  of  the  plant  model  and  compensator, 
respectively.  We  shall  assume1  that  the  plant  has  n  inputs  and  outputs,  and 
that  detP(s)  &  0.  The  signals  r(s),  d(«),  and  y(s)  are  the  reference  input, 
disturbance  input,  and  system  output,  respectively.  Define  the  input  open-loop 
transfer  function,  sensitivity  function,  and  complementary  sensitivity  function: 

LM  =  F(»)P(»),  s,(.)  =  [I  +  £,(»)|-\  JVM  =  i,(«)|J  +  £,(*)!"'  (3.2.1) 

and  the  output  open-loop  transfer  function,  sensitivity  function,  and  complemen¬ 
tary  sensitivity  function: 

L0[s)  -  P{s)F(s),  S0(s)  =  [I  +  L0{s)}-\  T0(s)  =  L0{s)[I  +  L0{s)}~1  (3.2.2) 

The  following  notation  will  be  used.  Let  RH°°  denote  the  set  of  proper 
rational  matrices  that  are  stable,  i.e.  that  have  no  poles  in  the  closed  right  half 
plane.  Given  G(s)  €  RH°°,  define  |jG||oo  =  supwff[G(;w)],  where  *[•]  denotes 
the  largest  singular  value  of  a  matrix  [GoV83].  Define  also  the  set 

D7  =  (<?(a)  6  Rff°°  :  IIGIU  <  7}  (3.2.3) 

When  we  refer  to  a  matrix  in  D7,  the  dimensions  will  be  clear  from  the  context. 
Finally,  we  suppress  dependence  upon  frequency  whenever  convenient. 

,Oor  remit*  extend  to  noneqaere  plant*  that  aatiefy  the  relevant  aeenmption  at  left  or  rifht 
invert  ibiHty. 
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Figure  3.1:  Feedback  Configuration 
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We  now  describe  the  design  specifications  that  our  system  is  to  satisfy.  Sup¬ 
pose  that  the  true  plant  is  given  by 

J*(.)  =  P(.)(/  +  fl,(.)A,(.)l  (SA.4) 

The  matrix  A;(i)  may  tie  anywhere  in  a  let  of  form 

A,  =  {  d*oy[Ai,  At] :  Ax  €  C*xk,  A,  €  A,  €  D7,  •  =  1,2} 

(3.2.5) 

where  0  <  k  <  n.  We  say  that  each  A<  is  an  unstructured  uncertainty.  The 
matrix  R/(s)  =  diag[ri(s)Ik,  rj(a)/n_*]  is  dimensioned  compatibly  with  A/(s) 
(/*  and  In-k  ve  identity  matrices).  The  weighting  functions  ri(a)  and  rj(a) 
are  chosen  to  be  stable  and  to  have  stable  inverses,  and  are  used  to  model  the 
frequency  dependence  of  the  level  of  uncertainty. 

Our  first  design  goals  axe  to  achieve  nominal  internal  stability  (cf.  [Vid85], 
Section  5.1)  and  to  maintain  stability  robust  against  the  uncertainty  described 
by  (3.2.4).  Using  standard  arguments  (e.g.,  see  [DoS81],  [DWS82]),  it  is  straight¬ 
forward  to  show  that  stability  will  be  robust  against  arbitrary  A<  6  D7  alone 
(i.e.,  with  Ay  =  0,  j  £  i)  if  and  only  if 

MflMl.  <  1/7  (3^.0) 

where  E\  =  H*  ,  Et  =  r°  .  Assessing  stability  robustness  against  both 

L  0  J  ln~h 

uncertainties  simultaneously  requires  use  of  the  structured  singular  value  [Doy82], 
and  will  be  deferred  until  Section  3.4. 

We  also  wish  to  satisfy  the  nominal  performance  goal  of  reducing  the  effect  of 
the  disturbance  input  upon  the  system  output.  We  say  that  this  goal  is  achieved 
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if 

\\r»S0\U  <  1  (3.2.7) 

where  r>(e)  is  a  weighting  function  whose  sise  indicates  the  relative  importance 
of  disturbance  reduction  at  each  frequency.  Next,  we  demand  that  performance 
be  robust  against  uncertainty.  Hence,  for  a  given  7,  we  require  that 

IM>|U<1,  VA/6A,  (3.2.8) 

where  S'Q  =  [I  +  PF]-1. 

We  shall  need  to  use  the  singular  value  decomposition  [GoV83]  of  the  plant  to 
describe  its  directionality  properties;  for  a  detailed  discussion,  see  e.g.,  [FrL86], 
[Fre80a],  and  [FVeOO].  Let  the  singular  value  decomposition  of  the  plant  trans¬ 
fer  function,  evaluated  at  a  fixed  frequency,  be  denoted  P  —  WTZHt  where 
T  =  diag[ tu  r,,  •••,  rj,  W  =  [u>i,  tea,  •••,»*],  and  Z  =  [su  *>,•••,*»].  The 
diagonal  elements  of  T  are  termed  the  singular  values  of  P,  and  are  ordered 
so  that  o[P]  =  fi  >  t%  >  -  •  ■  >  r*  =  g[P]  >  0.  The  columns  of  the  unitary 
matrices  W  and  Z  are  termed  the  left  and  right  singular  vectors,  respectively. 
The  largest  and  smallest  singular  values  of  P(j'w)  have  interpretations  as  the 
largest  and  smallest  possible  gains  at  frequency  u.  The  condition  number  of  the 
plant  is  defined  as  #c[P]  =  o\P\Iq[P\.  If  the  condition  number  is  very  large  at  a 
frequency,  then  the  gain  of  the  system  is  strongly  directional  at  that  frequency, 
and  the  plant  is  said  to  be  ill-conditioned  [GoV83j. 

It  will  prove  convenient  to  adopt  the  following  notation.  Consider  a  matrix 
whose  columns  form  an  orthonormal  basis  for  a  h-dimenskmal  subspace  of  C*. 
We  denote  the  matrix  by  an  uppercase  tetter  (or  a  lowercase  tetter  in  the  case 
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k  =  1),  and  the  associated  subspace  by  the  boldface  version  of  the  same  letter 
(e.g.,  the  columns  of  X  6  Cnxk  form  an  orthonormal  basis  for  the  lb-dimensional 
subspace  X  C  C“). 

Finally,  we  need  to  define  a  specific  class  of  compensators  which  we  shall  term 
plant-inverting  compensators.  Such  compensators  were  discussed  in  [fVe89b], 
[FreQO],  [SkM87],  and  are  useful  in  revealing  design  difficulties  associated  with 
ill-conditioned  plants.  A  plant-inverting  compensator  F(s)  is  of  the  form 

F(«)  =  l{s)P~l{s)  (3.2.9) 

where  l(s)  is  a  scalar  transfer  function  selected  so  that  F[s)  is  proper.  Obviously, 
use  of  a  plant-inverting  compensator  requires  that  the  plant  be  stable  and  have 
a  stable  inverse.  For  later  reference,  note  that  with  such  a  compensator  the 
two  sets  of  transfer  functions  (3.2.1-2)  are  identical.  We  shall  denote  these  by 
L0(s)  =  Lj(s)  =  /(«)/,  S0[s)  =  S/(«)  =  «(«)/,  and  T0{s)  =  Tt(s)  =  t(s)I,  where 
a(«)  a  1/(1  +  /(«))  and  t(«)  =  l(s)/(l  +  !(«)). 

3.3  First-  and  Second-Order  Effects  of  Input 
Uncertainty  upon  Output  Sensitivity 

The  robust  performance  problem  posed  in  the  preceding  section  requires  that  the 
output  sensitivity  function  be  kept  small  despite  the  presence  of  block  diagonal 
uncertainty  at  the  plant  input.  In  this  section  we  shall  study  this  problem  by 
investigating  the  first-  and  second-order  effects  of  input  uncertainty  upon  output 
sensitivity.  As  we  shall  show  in  Section  3.4,  an  understanding  of  only  the  first- 
and  second-order  effects  is  sufficient  to  allow  development  of  qualitative  design 
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rules  for  this  robust  performance  problem. 

Lemma  8.8.1  Assume  that  P  hat  the  farm  defined  in  (S.i.S-4)-  Then  S'0  = 
[/  +  PF]~l  hat  the  form  S'0  =  So  +  So  +  S$  +  O(Aj),  where 

=  -S0PRt*iP-lTo  (3.3.1) 

and 

Sq  =  SoPRibiTjRjbiP-'To  (3.3.2) 


Proof:  The  proof  follows  from  straightforward  manipulations  as  in  the  proof 
of  Lemma  4  of  [Re89a]. 

■ 

Lemma  3.3.1  allows  an  approximate  decomposition  of  S'0  into  the  sum  of  the 
seroth  order  term  So,  first  order  term  Sg,  and  second  order  term  Sg-  Achieving 
robust  performance  requires  that  Sg  and  Sg  be  kept  small.  Let  us  initially 
consider  the  first  order  effects. 

Proposition  8.8.3  Suppote  that  A /  6  A7.  Then 

rap  »[Sol  <l£>«l  ff[S0Pa]»[£fp-‘r0|  (J.s.3) 

A/€A*  <=1 

and 

sup  ^[5©]  >  7  max  |  j  9[SoPEi]o[E?  P~lTo]  (3.3.4) 

A/CA^  ^ 

Proof:  The  upper  bound  follows  by  noting  that  ■ 

So  =  —t\S0PE\&\E* P~lT0  -  rtSoPEiAtEfp-'Tg  (3.3.5) 
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and  applying  the  triangle  inequality.  The  lower  bound  follows  from  (3.3.5)  by 
noting  that 

■up  9[So]  >  sup  d(S^] 

Af€A«g  i  A|>0 

and  constructing  a  “worst  case”  Ai  as  in  the  proof  of  Lemma  5,  [Fre89a]. 


It  is  clear  that  the  terms*  |  r{  |  0[SoPEi\o[Ef  P~xTo\  characterise  the  first 
order  effects  of  uncertainty  at  the  plant  input  upon  the  output  sensitivity  func¬ 
tion.  Hence,  to  achieve  robust  performance,  it  is  necessary  that  keeping  these 
terms  small  be  made  a  design  objective. 

The  remainder  of  this  section  is  devoted  to  exploring  the  second  order  effects 
of  input  uncertainty  upon  the  output  sensitivity  function. 

Proposition  8.S.S  Suppose  that  A/  €  A,.  Then 

•up  <t’EE  I  nr,  1 9\SoPEl}B\E!'T,El)a\Bfp-1To)  (3.3.#) 

A/eA,  /=H*1 


Proof:  The  proof  follows  from  noting  that 

5*  =  r\S0PElAlEfTlE1^lEfP~lT0 

+  rxrtSoPEx^xEfTjE^tEfP^To  (3.3.7) 

+  rxrxSoPExbxEfTiEx&xE?P-xTo 
+  rlS0PEt&iEfTjEx&iEfP~xT0 

3 la  Section  3.4,  we  ahaD  allow  that  theee  Utrme,  u  well  ae  thooe  diecneeed  following  Propoehion 
3.3.8,  are  doeeljr  related  to  oar  interactum  prrn meter*. 
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and  applying  the  triangle  inequality. 


Proposition  3.3.3  shows  that  the  second  order  effects  are  characterised  fay  the 
terms  |  r<rj  1 9[SoPEt\9[EfTjEf\9[EfP~lTo]>  Hence  keeping  theee  terms  small 
should  also  be  included  as  a  design  objective.  In  this  regard,  we  notice  that 
the  first  and  fourth  terms  in  (3.3.7)  will  tend  to  be  small  if  the  robust  stability 
bounds  (3.2.6)  are  satisfied  and  if  the  first-order  effects  of  uncertainty  are  small. 
Analysis  of  the  second  and  third  terms  in  (3.3.7)  is  more  problematic,  and  we 
shall  return  to  this  point  at  the  cloee  of  Section  3.4. 

It  is  straightforward  to  extend  the  results  of  this  section  to  analyse  the  first- 
end  second-order  effects  of  input  uncertainty  modelled  as  an  mxm  block  diagonal 
matrix. 


3.4  Structured  Singular  Value  Analysis  of  the 
Robust  Performance  Problem 

The  first  step  in  applying  structured  singular  value  analysis  is  to  reformulate 
the  robust  performance  problem  as  an  equivalent  robust  stability  problem  by 
introducing  a  fictitious  modelling  uncertainty  to  represent  the  performance  spec¬ 
ification  [DWS82].  For  our  problem,  this  is  accomplished  by  supposing  that  the 
true  plant  is  given  by 

P(>)  =  (/  +  r,(.)A.(.)]-1J>(.)[/  +  *,(.)A,(.))  (3.4.1) 
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where  the  input  uncertainty  is  the  same  as  in  (S.2.3-5),  A*(a)  is  stable,  and  the 
weighting  function  rs(s)  is  identical  with  that  in  (3.2.7).  As  discussed  in  [Dcy82], 
[DWS82],  the  next  step  is  to  define 

A  =  diag[Au  As,  A,]  (3.4.2) 


and 


M  = 


'  rxE*TtEx 
rtE*TtEx 
—rtSoPEx 


rxExTiE%  —rxEfP~lT0 
r,  EfTjEt  -rtEfP-'To 
—r%SoPE%  ra5o 


(3.4.3) 


By  definition  [Dcy82],  the  structured  singular  value,  denoted  p[M(/(i;)],  is  the 


reciprocal  of  a  frequency  dependent  stability  margin  with  respect  to  all  three 


sources  of  uncertainty.  Define  jieo[AT]  =  supu  /i[M(jw)].  A  nominally  stable 


feedback  system  will  be  robustly  stable  against  all  uncertainties  of  the  form 
(3.4.2)  with  At  an  arbitrary  member  of  D,  if  and  only  if  #*oo[J \t]  <  1/7.  hi  terms 
of  the  original  robust  performance  problem,  we  have  that  lIrsSolloo  —  1/T  for  all 
A/  €  if  and  only  if  Hoo[M]  <  1  h  [DWS82].  For  the  rest  of  this  section,  we 
will  discuss  only  the  equivalent  robust  stability  problem. 


To  obtain  insight  into  the  robust  stability  problem,  we  will  work  with  a  set 
of  frequency  dependent  bounds  upon  the  structured  singular  value.  To  simplify 
notation,  we  denote  the  block  elements  of  M  by  and  define 


ay  =  o[Mij\ 


(3.4.4) 


and 


Pi  =  ma xou  (3.4.5) 

Hence  the  system  is  robustly  stable  against  the  uncertainty  A 4  alone  if  and  only 
if  <r«(ju>)  <  1/7,  Vw. 
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Next,  we  define  the  interaction  parameters  (cf.  Chapter  2) 


p(*.i)  =  y/Oi wit 

(3.4.6) 

(3.4.7) 

Ma  =  maxp(i,i) 

(3.4.8) 

P*  =  xnxxp{itj,k) 

(3.4.0) 

The  interaction  parameter  p(«,j)  determines  to  what  extent  the  effects  of  the 
two  uncertainties  A<  and  Ay  can  interact  to  produce  robustness  difficulties.  The 
robust  stability  problem  obtained  by  setting  A*  =  0  can  be  studied  using  struc¬ 
tured  singular  value  techniques  applied  to  the  matrix  obtained  by  deleting  the 
fcth  block  row  and  column  of  M.  Denote  this  matrix  by  M\ij.  The  following 
result  is  an  immediate  corollary  to  Proposition  2.3.1. 

Proposition  3.4.1  Choose  k.  Then 

max  {  on,  oyy,  #*(»,/) }  <  p[M\ij]  <  max  {  a#,  oyy  >  +  (3.4.10) 


For  our  problem,  the  interaction  parameters  in  (3.4.6)  are 

p(l,2)  =  V*MPZ>4]»Mf2V4j  (3.4.11) 

p(  1,3)  =  yJd\riSoPEi\o\r\EfP-lTo\  (3.4.12) 

p(2,3)  =  \lo[TtSoPEt)o[rtEfP-lTo]  (3.4.13) 

To  illustrate,  let  us  return  to  the  problem  of  robust  stability  against  block  di¬ 
agonal  uncertainty  at  the  plant  input.  Recall  that  the  system  is  robustly  stable 
against  A <  alone  (»  =  1,2)  if  and  only  if  (3.2.6)  is  satisfied.  The  system  will  be 
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robustly  stable  against  Ai  and  As  together  if  and  only  if  Moo[Af|i,j]  <  1/7*  It 
follows  from  Proposition  3.4.1  that,  in  addition  to  keeping  the  diagonal  blocks  of 
Tj  small,  it  is  also  necessary  to  keep  the  produet  of  the  off-diagonal  blocks  small. 
The  latter  requirement  is  to  prevent  destabilising  interactions  between  Ai  and 
As* 

The  interaction  parameters  ji(l,  3)  and  m(2,3)  each  have  an  analogous  inter¬ 
pretation.  Suppose  that  stability  is  robust  with  respect  to  Ai  and  As  individu¬ 
ally.  Then,  to  maintain  robust  stability  against  the  combined  effects  of  these  two 
uncertainties,  it  is  necessary  and  sufficient  that  m(1»3)  be  kept  small.  An  inspec¬ 
tion  of  Proposition  3.3.2  reveals  that  this  requirement  is  equivalent  to  keeping 
the  first  order  effects  of  Ai  upon  output  sensitivity  small.  Similar  comments 
apply  to  the  combined  effects  of  A*  and  As. 

We  next  consider  the  combined  effects  of  all  three  uncertainties.  The  following 
proposition  is  a  restatement  of  Propositions  2.3.2  and  2.3.3,  combined  here  for 
convenience. 

Proposition  3.4.2  Consider  Mi,  Ms*  and  ps  defined  in  (3. 4. 4-9).  Then 

(i) 

li[M)  >  max  {  mi,  ps,  Ms  >  (3.4.14) 

(ii)  Suppose  that  Ms  ^  0.  If  Ms/Ms  <  1»  Men 

fi[M\  <  Mi  +  M*  +  Ms  (3.4.15) 

(Hi)  Suppose  fit  ±  0.  If  Ms/ fit  <  1,  then 

M[M]  <  Mi  +  2.74ms  (3.4.16) 
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Far  later  reference,  we  note  that  for  our  problem, 


(*(1,2.3)  =  {?[riE?T,Et\t{rtSoPEl\»[TiEfp-'To§ V*  (S.4.1T) 

and 

(*(1, 3, 2)  =  (»(r,f:fr,f:1]»(r1S0PB,|»(rlf:fp-1r0|) 1/1  (3.4.18) 

From  an  inspection  of  Proposition  3.3.3  and  the  amming  discussion,  it  follows 
that  keeping  /ii,  p%  and  ns  ■«**!!  is  equivalent  to  keeping  the  serotb-,  first*,  and 
second-order  effects  of  input  uncertainty  upon  output  sensitivity  small.  Note  in 
particular  that  p(l,2,3)  and  p(l,3,2)  correspond  to  the  second  and  third  terms 
in  (3.3.7).  This  shows  that  the  essentials  of  the  robust  performance  problem  are 
captured  by  our  analysis  of  the  first-  and  second-order  effects  of  input  uncertainty 
upon  the  output  sensitivity  function. 

3.5  A  Robustness  Indicator 

In  the  previous  section  we  have  defined  interaction  parameters  that  must  be  kept 
small  if  performance  is  to  be  robust.  We  now  derive  a  robustness  indicator  whose 
sise  determines  the  potential  difficulty  of  this  task.  To  simplify  the  exposition, 
it  will  prove  convenient  to  assume  that  the  singular  value  decomposition  of  the 
plant  may  be  partitioned  as 

P  =  WXTXZ?  +  WxT,Zf  (3.3.1) 
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where  T\  =  d\ag[ru  r,,  •  •  • ,  ij],  T,  =  dtagfo+i,  Ti+J,  •  •  • ,  r„],  Wu  Zx  €  CBxl,  and 
Ws,  ^  €  Cnx(n-^.  We  shall  further  assume  that  rj  >  rj+ll  and  refer  to  the  plant 
as  consisting  of  “high  and  low  gain  subsystems”  [FrL86]  with  input  directions 
Zj  and  output  directions  W|.  Finally,  we  shall  assume  that  the  gain  in  each 
subsystem  is  uniform:  n  =  q,  and  tj+i  =  rn.  If  any  of  the  above  assumptions  fail 
to  hold,  analysis  may  be  performed  using  the  dyadic  form  of  the  plant  singular 
value  decomposition,  and  steps  toward  accomplishing  this  will  be  indicated  in 
the  sequel. 


We  next  need  to  describe  how  well  the  plant  input  directions  are  aligned  with 
the  two  blocks  of  uncertainty.  We  do  tins  by  introducing  an  angular  measure 
of  the  distance  between  the  subspaces  Zi  containing  the  inputs  to  the  high-  and 
low-gain  subsystems  of  the  plant  and  the  subspaces  E|  containing  the  inputs  to 
and  outputs  from  each  block  of  uncertainty.  Consider  the  subspaces  Zi,  Zi,  Ei, 
and  £],  and  let  q  =  min{k,  /}.  As  discussed  in  Section  2  of  [FrL86]  and  the 
references  cited  therein,  the  principal  angles  [BjG73],  [GoV83]  &  =  otj  >  as  > 
•  •  ■  >  a,  =  e,  a*  G  [0,  jt/2],  between  the  subepaces  Zi  and  Ei  are  a  measure  of 
the  angular  distance  between  these  two  subspaces,  and  thus  also  between  their 
orthogonal  complements  Z)  and  Ea.  These  angles  have  the  following  useful 
characterisations. 

Lemma  8.S.1  Consider  the  matrix 


E*Z  = 


EfZ\  EfZt 
EfZx  Ef  Zx 


(3.5.2) 


Assume8  that  k+l  <  n.  Then  there  exist  compatibly  dimensioned  unitary  matrices 

9Otherwise  (n  -  1)  +  (n  -  k)  <  n,  and  the  decomposition  may  be  obtained  by  renumbering 
and  replacing  k  with  n  -  k  and  l  with  n  —  l. 
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U  =  diag[Uu  Ut]  andV  =  diag[Vu  Vt]  such  that 
( i)  if  k>  l,  then 


U*E?ZV  = 


C  :  -S  0 

0  :  0  /*_, 


0 

0 


5  :  C  0  0 

.  0  :  0  0  /*-(**) 


(3^.3) 


(it)  if  k  <  l,  then 


C  0  :  -S  0 


UB1 &*ZV  = 


S  0 
0  h-k 
0  0 


C  0 

0  0 

0  i»-(k+ 1) 


(3.5.4) 


where  C  —  diag\<umau  co 8aSl  •••,  cos  a,]  and  S  s  diag[amau  ainaj,  •••,  ainat,]. 


Proof:  A  straightforward  generalisation  of  the  C~S  Decomposition  of  a  general 
unitary  matrix,  which  is  presented  in  [Ste77],  [GoV83]. 

(Comparing  (3.5.2)  with  (S.5.3-4),  it  follows  that  the  singular  values  of  EfZy 
are  equal  to  the  cosines  of  the  principal  angles,  while  the  singular  values  of  EfZ% 
are  equal  to  the  sines  of  the  principal  angles  and  (possibly)  unity,  and  so  forth. 
This  proves  significant,  as  it  follows  from  (3.3.1)  and  (3.5.1)  that  the  first-order 
effects  of  input  uncertainty  upon  the  output  sensitivity  function  depend  upon 
both  the  plant  condition  number  and  the  matrices  EfZj,  i,  j  =  1,3.  The 
following  lemma  makes  this  statement  precise. 

Lemma  8.5.3  Suppose  that  the  plant  is  partitioned  as  in  ($.5.1).  Then,  from 
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(8.8.1),  it  follow  that 


Sb  =  -rxSoW&ZfEtAiEfZiTT'WfTo 

-rtSoW&ZfEtAiEfZtTZ'WfTo  (3.5.5) 

+other  terms 

Note  that  the  first  two  terms  in  this  expansion  are  proportional  in  sise  to  the 
plant  condition  number.  Hence,  if  the  plant  is  ill-conditioned,  the  first-order  ef¬ 
fects  in  (3.3.1)  are  potentially  large.  By  the  uniform  gain  assumption,  the  “other 
terms1'  in  this  expansion4  are  not  proportional  to  the  plant  condition  number, 
and  thus  are  of  less  significance.  Also,  this  result  shows  that  the  extent  to  which 
plant  ill-conditionedness  can  cause  design  difficulty  depends  upon  the  alignment 
between  the  subspaces  Ef  and  Zj.  To  illustrate,  consider  Sq  as  the  series  con¬ 
nection  of  S0P ,  .R/A/  and  P~lT0  (c.f.  (3.3.1)).  Potential  design  difficulties  arise 
whenever  the  uncertainty  can  introduce  coupling  between  the  high-  (km-)  gain 
subsystems  of  P  and  P~l.  For  example,  if  k  =  l  =  n/2,  then  no  coupling  will  be 
introduced  if  the  alignment  is  perfect  (i.e.,  ft  =  0  or  g  =  w/2).  Certain  combir 
nations  of  the  dimensions  k,  l,  and  n  will  always  result  in  coupling,  regardless  of 
how  the  subspaces  are  aligned.  For  an  example,  consider  k  =  1,  l  —  2,  and  n  =  4. 
Also,  a  limiting  case  is  when  k  =  n ,  for  which  the  uncertainty  is  no  longer  block 
diagonal  and  the  analysis  reduces  to  that  of  [Fre89aj.  The  potential  robustness 
difficulties  arising  from  these  and  other  combinations  of  k,  l ,  and  n  can  be  best 
summarised  by  considering  the  interaction  parameters  /i(l,3)  and  p{ 2,3). 
Proposition  3.5.3  Suppose  that  the  partition  (8.5.1)  and  the  uniform  gain  as - 

4It  it  straightforward  to  calculate  these  terms  explicitly.  Since  the  resulting  expressions  are 
lengthy,  they  are  omitted  for  purpose  of  brevity. 
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sumption  hold,  and  that  k  +  l<n.  Then 


m’(I.S)  >  e[SoWroMi>l  (3-6.®) 

<**(a,s)  >  slSoJeiroWli’)  (sat) 

where 

t<! 

sin &,  k  +  l  <n,  k>l 

_ rpi  _  .  |  r  pi  CM  A  Bin  fi,  k  +  l  =  n,  k>l 

Pu[P)  -  |  r*r$  |  *[P]  lf  *  +  /<n|jk</ 

cm g,  k  +  l  =  n ,  k<l 

Furthermore,  for  i  =  1, 2, 

M*(«.  3)  <  r0]A»[P]  +  other  terms  (3.5.10) 


Proof:  Using  standard  norm  inequalities  and  the  fact  that  W  is  a  unitary 
matrix  yields 

(*’(1.3)  >1  r,r,  |  <dS0\e\T0)»{T<Z!'E1ME?ZlTfl\ 

Setting  *  =  1,  j  =  2  and  applying  the  uniform  gain  assumption  gives 

M*(l,3)  >|  r,r,  |  <i[So\<i{T0\k[P\b\Z* Ev\t\E? Z,\, 

from  which  (3.5.6)  follows.  Proof  of  (3.5.7)  is  similar.  The  upper  bound  (3.5.10) 
follows  fay  substituting  (3.5.1)  into  (3.4.12-13)  and  applying  norm  properties. 


(3.5.8) 

(3.5.0) 
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It  follows  from  (S.5.6-7)  that  robustness  Is  potentially  poor  at  any  frequency 
for  which  the  robustness  indicator 

p[P]  =  maxp<a[P]  (3.5.11) 

is  large.  At  such  frequencies,  it  is  necessary  that  the  output  loop-transfer  func¬ 
tion  be  shaped  so  that  p[5n|o[2V)]  is  small.  Since  the  "other  terms*  in  (3.5.10) 
are  not  proportional  to  the  plant  condition  number,  it  follows  that  #*(*,  3)  may  be 
kept  ntnii.11  by  requiring  the  product  0\SoW^O^/^To\  to  be  small.  Comparison 
with  (3.5.5)  shows  that  this  strategy  is  equivalent  to  keeping  the  first-order  ef¬ 
fects  of  Ai  and  Aj  small.  For  purpose  of  comparison  with  the  results  of  [Fke80a], 
[Fre89b],  [FraOO],  and  [Ste87]  we  note  that  the  case  of  unstructured  input  uncer¬ 
tainty  corresponds  to  setting  k  —  n  (3.2.5).  In  this  case  the  weighting  rs  is 
irrelevant  and  6,  =  0,  and  hence  the  robustness  indicator  assumes  a  correspond¬ 
ingly  simple  form: 

p[P]  =|  rir,  |  k[P\  (3.5 J 3) 

which  depends  solely  upon  the  plant  condition  number. 

In  summary,  we  have  shown  that  if  the  robustness  indicator  (3.5.11)  is  large, 
then  performance  robustness  is  potentially  poor.  The  analysis  is  limited,  how¬ 
ever,  in  that  it  is  based  solely  upon  the  interaction  parameters  /t(i,3),  “d  thus 
solely  upon  the  first  order  effects  of  uncertainty.  It  remains  to  be  determined 
whether  performance  robustness  can  be  potentially  poor  because  of  the  second- 
order  effects  of  uncertainty,  even  if  the  first-order  effects  are  negligible.  This 
question  translates  approximately  into  that  of  determining  whether  the  interac¬ 
tion  parameters  p(i,j,fc)  can  be  large  even  if  H\  and  are  small.  A  somewhat 
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surprising  result,  which  is  applicable  to  systems  with  two  inputs  and  two  outputs, 
shows  that  keeping  the  interaction  parameters  and  /* %  small  tends  to  insure 
that  /is  is  small  as  well. 

Proposition  3.5.4  Let  n  =  2,  and  let  fii,  /it  and  Ms  be  defined  in  (S.4.4-9). 
Then 

Ms  <  1.62  max  {  Mt,  Mi }  (3.5.13) 

Proof:  See  Appendix  A. 

■ 

Suppose  that  the  goals  of  nominal  performance  and  robust  stability  against 
diagonal  input  uncertainty  ere  satisfied,  so  that  Mi  is  small.  Then  it  follows 
from  Propositions  3.4.2  and  3.5.4  that  the  goal  of  achieving  robust  performance 
essentially  translates  into  that  of  keeping  Ms  small.  It  follows,  therefore,  that  the 
robustness  indicator  does  adequately  reflect  the  potential  difficulty  of  achieving 
robust  performance.  Finally,  although  it  seems  reasonable  to  conjecture  that 
Proposition  3.5.4  can  be  extended  to  the  general  case  of  block  2x2  uncertainty, 
the  proof  of  this  proposition  does  not  readily  generalise,  and  the  conjecture  is 
left  for  future  research. 

In  what  follows  we  shall  discuss  the  robustness  difficulties  associated  with  a 
plant-inverting  compensator.  Note  first  that  the  interaction  parameter  Ms  =  0  in 
this  case,  so  that  the  discussions  following  Proposition  3.5.4  are  rendered  moot. 
To  explain  these  difficulties,  we  first  recall  the  following  result  of  Stein  [Ste87], 
who  explicitly  calculates  the  structured  singular  value  resulting  from  use  of  a 
plant-inverting  compensator  in  the  case  of  unstructured  uncertainty. 
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Proposition  8.5.5  [Ste87]  Assume  that  k  =  n  in  (8.8.5),  and  that  F(a)  is 
given  by  (8.8.0).  Then 

H*[M]  -  M*  +  |r,a|*  +  Ifitj  •  \rts\  •  (*[P]  +  l/ic[P])  (3J5.14) 


Hence,  it  follows  that,  for  a  fixed  loop  shape  /(«),  the  structured  singular  value 
will  increase  as  the  square  root  of  the  robustness  indicator  (3.5.12). 


We  next  develop  a  corresponding  result  for  the  present  problem.  For  no- 
tational  simplicity,  we  consider  only  the  special  case  k  —  l  =  n/2,  for  which 
p[P\  =  max*  |r<rs|/c[P]  coefiainfi.  Noticing  that  p s  =  0,  /*(1,2)  =  0,  the  bounds 
given  in  (3.4.14)  and  (3.4.16)  can  be  utilised.  However,  stronger  results  may  be 
obtained  due  to  the  special  structure  of  the  interconnection  matrix. 
Proposition  8.5.6  Assume  that  k  =  l  =  n/2  in  (8.8.5),  and  that  the  partition 
(8.5.1)  and  the  uniform  gain  assumption  hold.  Let  F(a)  be  given  as  in  (8.8.0). 
Then 

max  {pi,  mi}  <  m(M]  <  Mi  +  V^Ma  (3.5.15) 


where 

Mi  =  max  {|rxt|,  |r,t|,  |r,«|} 

Mi  =  max  {m(1,3),  m(2,3)> 

Furthermore, 

M*(»,  3)  >  |r<t|  *  |rs«|ic[P]cosflsinfi 


(3.5.16) 


(3.5.17) 


and 


m’(*«3)  <  Ml  *  |r*«|  (a(P]cosfisinA  +  coeafl  +  sin,a+  (l/a(P])cosflsintt) 

(3.5.18) 
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Proof:  See  Appendix  B. 


Hence  again,  we  see  that  for  a  fixed  loop  shape,  the  structured  singular  value 
will  increase  with  the  square  root  of  the  robustness  indicator  (3.5.11).  The  main 
difference  is  that  now  the  robustness  indicator  depends  not  only  upon  the  plant 
condition  number,  but  also  upon  the  alignment  angles.  Motivated  by  these  anal¬ 
yses  and  those  in  [SkM87],  it  is  of  particular  interest  to  develop  an  intuitive 
explanation  of  the  design  difficulties  associated  with  such  a  compensator.  To¬ 
ward  this  end,  note  that  when  the  plant  has  the  form  (3.5.1),  such  a  compensator 
also  has  a  decomposition  into  high  and  low  gain  subsystems.  In  the  absence  of 
uncertainty,  the  high  (or  low)  gains  of  the  compensator  cancel  the  low  (or  high) 
guns  of  the  plant  exactly.  However,  with  uncertainty  present,  this  cancellation 
will  be  exact  only  i/the  alignment  conditions  are  satisfied  (and  hence  the  robust¬ 
ness  indicators  are  sero).  Otherwise,  the  uncertainty  can  cause  the  high  (or  low) 
guns  of  the  plant  and  compensator  to  multiply  one  another,  so  that  the  resulting 
series  connection  of  plant,  uncertainty,  and  compensator  is  significantly  different 
from  its  nominal  value.  Remarks  similar  to  theae  are  presented  in  [SMD88],  and 
physical  interpretations  are  presented  for  robustness  problems  associated  with 
high-purity  distillation  columns.  Our  results  above  are  appealing  in  that  they 
allow  these  potential  difficulties  to  be  quantified  in  terms  of  the  principal  angles. 

Finally,  if  the  uniform  gain  assumption  fails  to  hold,  then  analysis  can  pro¬ 
ceed  using  the  dyadic  form  of  the  plant  singular  value  decomposition:  P  — 
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r*w**F-  Indeed,  calculations  similar  to  those  in  Proposition  3.5.3  yield 

/**(•’. 3)  >|  r<r»  |  s[So]s[To\  •  {Tj/Tk)\\zf  Ei\\  •  ||E<"  **||  (3.5.19) 

where  ||sjFE<||  and  ||.E<ir**||  are  the  cosines  of  the  angles  between  the  subspaces 
Sj  and  £i,  and  s^  and  E{,  respectively.  Hence  p[P]  =  max<pa[P],  where 

*»[P]  =  I  rtr»  |  (ty/rfc)||*f  X||  •  ||^s*||  (3.5.20) 

is  a  useful  generalisation  of  the  robustness  indicator  to  this  more  general  situa¬ 
tion. 

3.6  The  Relative  Gain  Array 

The  Relative  Gain  Array  (RGA)  [Bri66]  is  widely  used  in  the  process  control 
industry  to  analyse  interactions  in  multivariable  systems.  A  number  of  authors 
have  studied  robustness  difficulties  associated  with  plants  having  large  elements 
in  the  RGA.  (For  discussions  of  these  results  and  lists  of  references,  see  [NeM87], 
[SkM87].)  Skogestad  and  Morari  [SkM87]  have  demonstrated  that  the  RGA  is 
useful  as  an  indicator  of  potential  design  difficulty  for  the  problem  of  maintain¬ 
ing  performance  robust  against  diagonal  multiplicative  uncertainty  at  the  plant 
input.  In  particular,  they  show  that  the  RGA  is  a  less  conservative  indicator 
than  is  the  plant  condition  number  for  this  particular  problem.  On  the  other 
hand,  they  also  show  that  the  RGA  can  be  optimistic,  and  can  sometimes  fail 
to  detect  potential  difficulties  when  the  plant  has  large  off-diagonal  elements  (cf 
[SkM87]).  Our  purpose  in  this  section  is  to  compare  our  robustness  indicator 
with  the  RGA,  and  to  show  that  ours  does  not  suffer  from  thin  shortcoming.  It 
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suffices  to  consider  the  case  n  =  2,  and  to  assume  that  the  weightings  in  (3.2.4) 
are  unity. 


For  our  purpose,  the  SGA  is  a  matrix 

RGA  = 


An  A« 

An  Am 


(3.6.1) 


where  A»  =  An,  AM  =  An  =  1  -  An,  Aa  =  1/  (l  -  *“£)  and  p ij  =  1,2  are 
the  elements  of  the  plant  transfer  function  matrix.  For  n  =  2,  the  subspaces  Zf 
and  E(  are  one-dimensional  and  thus  will  be  denoted  as  i|  and  ei,  here  and 


correspond  to  the  vectors  — 


and  E)  = 


o' 

1  J‘ 


p[P]  =  /c[P]  cos  a  sin  a 


The  robustness  indicator  is 
(3.6.2) 


here  a  denotes  the  principal  angle  between  the  subspaces  si  and  ei.  The  following 
is  adapted  from  [SkM87|. 

Proposition  3.6.1  Suppose  that  uncertainty  has  the  form  of  (9.8.1)  with  n  = 
2,  k  =  1.  Then  the  output  open-loop  transfer  function  eon  be  written  as  VQ  = 
(I  +  E)PF,  where  E  is  an  error  term: 

E  =  PA/P"1  (3.6.3) 


Furthermore, 


sup  9[E]  >  7  max  I  Au  |*|  {—  | 
A/e a,  *J  Pa 


Proof:  The  proof  follows  from  noting  that  (SkM87j 


PA/P-1  = 


AuAi  +  AjjAj  -Au(pit/pn)(Ai  -  As) 
[  -Au(pn/pii)(Ai  -  A«)  AuAi  +  Am  As 


(3.6.4) 


(3.6.5) 
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It  follows  from  (3.6.4)  with  *  =  j  that  large  values  of  the  relative  gain  |  An  | 
correspond  to  systems  that  are  very  sensitive  to  uncertainty.  However,  small 
values  of  |  An  |  need  not  necessarily  correspond  to  systems  that  are  insensitive. 
Indeed,  counterexamples  can  readily  be  constructed  using  a  triangular  P  matrix 
whose  largest  element  is  on  the  off-diagonal  (e.g.,  p.  2328,  [SkM87]). 

Expressing  the  plant  output  with  respect  to  the  basis  consisting  of  the  left 
singular  vectors  of  the  plant  displays  the  role  played  by  the  condition  number 
and  input  directionality  properties  of  the  plant. 

Proposition  S.6.2  Consider  the  error  term  (8.6.8): 

sup  8[E]  >  'imax{Ti/rk)  |  sfe,  |  •  |  efs*  |  (3.6.6) 

a,eA, 

and 

sup  0\E]  <  4Tfm«(t)/r*)  |  *f  e<  |  ■  |  ef*k  |  (3.6.7) 

A  ,€A, 


Proof:  The  proof  follows  by  applying  norm  properties  to  the  error  term  J[E], 
and  noting  that  8{E\  =  8tyVBEW\%  where 


WbEW  = 


E?=i  si  KfPjEfsi^A^efsi 
(l/^l^])  E?=i  z?ei&ic?zl  2j.i  *fe<  A<ef  x. 


(3.6.8) 


In  particular,  it  follows  from  (3.6.2)  and  (3.6.6)  that 

sup  fr[E]  >  lp[P]  (3.6.9) 

A,eA, 
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Finally,  the  upper  bound  (3.6.7)  shows  that  the  error  term  cannot  be  significantly 
larger  than  0[A/]  if  the  robustness  indicator  is  small 


3.7  Examples 


The  authors  of  [SMD88]  illustrated  the  design  difficulties  of  ill-conditioned  plants 
using  simplified  models  of  the  LV  and  DV  configurations  of  a  distillation  column. 
Both  configurations  result  in  an  ill-conditioned  plant,  and  both  display  robustness 
problems  when  a  plant-inverting  compensator  is  used  and  the  input  uncertainty 
is  assumed  to  be  unstructured.  If,  however,  the  input  uncertainty  is  assumed  to 
be  diagonal,  then  only  the  LV-configuration  continues  to  pose  robustness  prob¬ 
lems  when  a  plant-inverting  compensator  is  used.  Based  upon  an  analysis  of 
the  directionality  properties  of  the  plants,  the  authors  of  [SMD88]  presented  an 
intuitively  appealing  explanation  of  this  phenomenon.  Our  intent  in  this  section 
is  to  demonstrate,  using  our  robustness  indicator,  that  the  DV-configuration  is 
indeed  less  difficult  to  robustly  control  when  the  input  uncertainty  is  diagonal. 


The  plant  models  for  both  configurations  are  given  by 


Plv(?)  ~  (1/(75*  +  1)] 


0.878  -0.864 
1.082  -1.096 


(3.7.1) 


and 


Pnv  =  [l/(75s  +  l)] 


-0.878  0.014 
-1.082  -0.014 


(3.7.2) 


respectively.  The  condition  numbers  are  constant  with  frequency:  k[Plv]  = 
141.7,  k\Pdv)  =  70.8.  The  weighting  functions  used  in  [SMD88]  are  rt(s)  = 
r*(a)  =  r(«)  =  0.2(5*  + 1)/(0.5«  +  1),  and  r*(«)  =  0.5(10*  +  l)/10a.  These 
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data  readily  verify  that  the  robustness  indicator  for  the  case  of  unstructured 
uncertainty  (cf  (3.5.12))  is  large.  By  Proposition  3.5.5,  a  plant-inverting  com¬ 
pensator  will  therefore  yield  poor  robustness  properties,  and  this  was  indeed 
verified  in  [SMD88]  using  Bode  plots  of  the  structured  singular  value.  When  a 
diagonal  uncertainty  structure  is  assumed,  however,  the  robustness  indicator  is 
given  by  (3.6.2).  The  principal  angles  for  the  two  configurations  are  a^v  =  45° 
and  (toy  =  0°,  respectively.  Hence  only  the  LV-configuration  poses  potential  ro¬ 
bustness  difficulties,  a  fact  also  verified  in  [SMD88]  using  the  structured  singular 
value  plots. 


To  further  illustrate  the  dependence  of  the  structured  singular  value  upon 
the  plant  input  directions,  we  alter  the  directions  of  the  LV-configuration  and  in 
the  meanwhile  keep  its  condition  number  unchanged.  Define 


P#(«)  =  Plv(*)R, 


R  —  Z 


cos  4  —sin  4 
sin  4  cos  4 


(3.7.3) 


where  Z  = 


cos  45°  -sin  45° 
sin  45°  cos  46° 


is  a  matrix  of  the  right  singular  vectors  for 


PlvUu)  [SMD88].  Hence,  the  robustness  indicator  for  this  new  plant  is 


p[Pt\  —  |rr8  |«[PlvJ  cos  4  sin  4  (3.7.4) 

Plots  of  the  structured  singular  value  for  the  compensator*  F(s)  =  (0.7 /s)Pt~1(e) 
versus  different  values  of  4  are  shown  in  Figure  3.2,  and  indeed  we  see  that  these 
plots  are  consistent  with  the  analyses  in  Section  3.5. 

Finally,  we  note  that  the  problem  of  aligning  the  plant  input  directions  is  one 

•Note  that  this  compensator  was  discussed  in  [SMDMj  for  6  =  45°  and  when  the  input 
uncertainty  is  unstructured. 
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magnitude  (absolute) 


frequency  (rad/sec) 


Figure  3.2:  SSV  Plots  for  Plant-Inverting  Compensators  with  Respect  to  P§. 
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of  plant  configuration  rather  than  that  of  compensation.  This  follows  because  the 
Input  uncertainty  is  assumed  to  be  located  between  the  plant  and  compensator. 

3.8  Concluding  Remarks 

In  this  Chapter  we  have  studied  a  robust  performance  problem  that  arises  when 
the  transfer  function  of  plant  model  is  ill-conditioned  and  the  mntfallitig  uncer¬ 
tainty  at  the  plant  input  has  block  diagonal  structure.  Our  analysis  shows  that 
design  goals  should  include  keeping  the  first-  and  second-order  effects  of  input 
uncertainty  upon  the  output  sensitivity  function  small  in  order  to  achieve  robust 
performance,  and  this  poses  potential  design  difficulties  for  plants  having  large 
robustness  indicators. 
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Chapter  4 

Several  Sources  of  Modelling 
Uncertainty 


4.1  Introduction 

An  important  problem  in  feedback  design  is  to  insure  robustness  of  stability  and 
performance  with  respect  to  several  sources  of  modelling  uncertainty.  A  useful 
stability  margin  for  an  important  class  of  modelling  uncertainties  can  be  cal¬ 
culated  using  Doyle’s  structured  singular  value  (SSV)  [Doy82],  [DWS82].  By 
plotting  the  structured  singular  value  as  a  function  of  frequency,  robustness  diffi¬ 
culties  associated  with  a  given  design  can  be  reliably  detected.  However,  because 
the  structured  singular  value  is  a  complicated  function  of  the  plant,  compensator 
and  uncertainty  description,  it  is  sometimes  difficult  to  obtain  design  insights 
using  only  such  a  plot.  For  example,  if  the  system  has  a  poor  stability  margin, 
one  would  like  to  know  how  this  may  be  improved  by  adjusting  the  compensator. 
Moreover,  some  feedback  design  problems  are  inherently  difficult,  and  one  would 
also  like  to  know  whether  one  is  faced  with  such  a  problem  and,  if  so,  how  a 
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compensator  may  be  selected  to  achieve  a  judicious  tradeoff  among  conflicting 
design  objectives.  The  computation  of  the  structured  singular  value  via  numer¬ 
ical  optimisation  (cf.  [FaT86],  [FaT88],  [Hel88],  [PFD88])  renders  such  insights 
difficult  to  obtain. 

Our  goal  in  the  present  paper  is  to  develop  results  that  should  be  useful  in 
obtaining  insights  into  robustness  problems  associated  with  several  modelling 
uncertainties.  To  motivate  our  approach,  let  us  consider  a  procedure  whereby 
we  assess  stability  robustness  sequentially,  first  with  respect  to  uncertainties 
taken  one  at  a  time,  then  taken  two  at  a  time,  and  so  forth.  As  we  proceed, 
we  may  find  that  a  system  which  is  robust  against  all  combinations  of  Jb  or 
fewer  uncertainties  can  nevertheless  be  destabilised  by  combinations  of  k  +  1, 
possibly  small,  uncertainties.  Evidently,  such  an  instability  must  be  caused  by  a 
harmful  interaction  among  these  k  + 1  uncertainties.  It  is  therefore  of  practical 
interest  to  have  a  methodology  for  determining  the  existence  and  strength  of  such 
interactions.  Our  analysis  below  will  identify  a  set  of  interaction  parameters, 
which  appear  to  be  useful  for  this  task. 

Essentially,  our  approach  in  this  chapter  is  a  generalisation  of  that  taken  in 
Chapters  2  and  3.  We  provide  additional  interpretations  of  the  earlier  results,  and 
lay  the  groundwork  for  applying  our  general  approach  to  other  design  problems. 

The  remainder  of  this  chapter  is  outlined  as  follows.  In  Section  4.2,  we  review 
the  definition  and  properties  of  the  structured  singular  value.  In  Section  4.3  we 
define  the  interaction  parameters  and  show  they  are  related  to  Mason’s  gun  rule. 
In  Section  4.4  we  derive  bounds  upon  the  structured  singular  value  expressed  in 
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terms  of  these  parameters.  Section  5  contains  an  analysis  of  the  aforementioned 
robust  performance  problem,  and  we  conclude  our  discussion  in  Section  4.6.  Some 
of  the  results  described  in  this  chapter  are  discussed  in  [ChF89b]  and  [ChFOO]. 


4.2  Preliminaries:  Properties  of  the  Structured 
Singular  Value 


The  first  step  in  the  structured  singular  value  analysis  is  to  rearrange  the  block 
diagram  of  the  feedback  system  into  the  form  shown  in  Figure  2.1.  The  uncer¬ 
tainty  matrix  A  (a)  is  assumed  to  lie  in  the  set1 


A.,={4iiay[A<(*)] :  A <  €  C"!*"*  and  stable,  || A^lfoo  <  7,  »  =  1,2,  •  •  •  ,N} 

(4-2.1) 

where  ||  Ai||«,  =  supw  ff[A<(j'w)],  and  ?[•]  denotes  the  largest  singular  value,  or  the 
Euclidean  norm  of  a  matrix.  The  interconnection  matrix  M[»)  is  a  function  of 
the  nominal  plant  model,  the  compensator,  and  weighting  matrices  introduced 
to  normalise  the  size  of  the  uncertainties.  It  is  stable  if  the  feedback  system 
is  nominally  stable.  At  each  frequency,  the  structured  singular  value,  denoted 
/i[M],  is  defined3  by  [Doy82]. 


0 

l/(minA€A..{0[A] :  det[I  +  M Aj  =  0}) 


det[I  +  M Aj^OVAG  A« 
otherwise 

(4.2.2) 


1  We  imuiu  throughout  this  chapter  that  all  A/s  in  (4.2.1)  an  distinct.  Tha  general  case  at 
repeated  A/s  ia  treated  in  [Doy82],  [P7DM|. 

aWe  shall  suppress  dependence  upon  frequency  whenever  appropriate,  and  shall  occasionally 
abuse  notations  by  using  the  symbol  A7  to  denote  the  set  of  constant  complex  matrices  with 
block  diagonal  structure  (4.2.1). 
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Define  #*oo[M]  =  sup*  fi[M(ju)}.  Then,  if  the  feedback  system  in  Figure  2.1  is 
nominally  stable,  it  will  be  robustly  stable  against  all  uncertainty  in  the  set 
if  and  only  if  [M]  <  1/7. 

Further  properties  of  the  structured  singular  value  are  summarised  in  the 
following;  for  a  complete  discussion,  see  [Doy82].  Define  the  sets 

D  =  {  diag[dxIni ,  •  •  • ,  <**/„„]  ‘  d*  €  R,  di  >  0  }  (4.2.3) 

where  denotes  the  identity  matrix  of  dimension  n*,  and 

U  *  {  diag[ «i,  U),  •  •  • ,  it*] :  ||u<||  =  1,  «<  €  C"»xl  >  (4.2.4) 

The  following  facts  provide  tools  for  evaluating  the  structured  singular  value. 
Fact  4.2.1  [Doy82] 

p[M\  <  b ri9[DMD~l)  (4.2.5) 

Fact  4.2.2  [Day82]  Let  p[-|  denote  the  spectral  radius  of  a  matrix.  Then, 

u^t{XaMY\  (4A8) 

Note  in  particular,  the  upper  bound '  (4.2.5)  is  equal  to  p[AT]  when  the  system 
has  three  or  fewer  blocks  of  uncertainty,  namely,  N  <  3. 

We  shall  partition  the  interconnection  matrix  into  blocks  whose  dimensions 
are  compatible  with  those  of  the  uncertainty: 

M(.)  =  [«,(.)],  i,jrn  1, N,  €  C**^  («.J) 

and  define  =  0[hUj}.  It  will  often  be  convenient  to  consider  the  simplified  sta¬ 
bility  robustness  problem  obtained  by  ffl^iiming  that  a  subset  of  the  uncertainty 


50 


blocks  A i  are  equal  to  sero.  We  introduce  the  following  notation  to  describe  such 
a  problem.  Let  Zn  =  { 1,  2,  •  •  • ,  N  },  and  define 


N*  =  {{*i»  *  *  *  *  •*}  ’  *i»  •*•»•*€  Zji  ,  «i  £  •  •  •  /  **}  (4.2.8) 

The  cardinality  of  K.  u  (  ^  W.  ^  the  <*  K* 

have  been  ordered  in  some  fashion,  and  shall  denote  these  elements  by  Jw, 
t  =  1,2,* ^  ^  Hence  for  each  »,  Jw  is  a  subset  of  Zu  with  k  distinct 
elements  (see  Example  4.3.1  for  an  illustration).  Define 

=  {A  :  A  €  A,,  A,  =  0,  V/  *  J«}  (4.2.9) 


and 


M\ju  =  {interconnection  matrix  M  (4.2.7)  with  the  jtk  block  row 

and  column  set  equal  to  sero,  Vj  £  Jw.}  (4.2.10) 

Then,  it  follows  readily  that  the  feedback  system  pictured  in  Figure  2.1  is  robustly 
stable  against  all  uncertainties  of  the  form  (4.2.9)  if  and  only  if  Moo[M|/m]  <  1/7. 
For  example,  suppose  that  only  one  of  the  A/s  in  (4.2.1)  is  nonsero;  i.e.,  suppose 
that  we  consider  the  set  Ju  =  {*'}.  Then,  it  follows  from  the  definition  (4.2.2) 
that  the  feedback  system  in  Figure  2.1  is  robustly  stable  against  all  uncertainties 
in  the  set  A,|/M  if  and  only  if  o«(jw)  <  1/7,  Vw. 

4.3  Interaction  Parameters 

From  the  analysis  at  the  close  of  the  preceding  section,  it  follows  that  if  the 
only  nonsero  blocks  of  the  interaction  matrix  are  those  on  the  diagonal,  then 
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p[M]  —  m *x<  an.  In  general,  however,  nonsero  off-diagonal  blocks  will  introduce 
interactions  between  blocks  of  uncertainty,  and  thus  will  affect  the  else  of  the 
structured  singular  value.  In  this  section  we  introduce  a  set  of  interaction  pa¬ 
rameters  whose  magnitudes  are  a  function  of  the  off-diagonal  blocks  of  M,  and 
we  shall  argue  that  these  parameters  determine  the  extent  to  which  interactions 
among  two  or  more  blocks  of  uncertainty  can  cause  stability  robustness  problems. 

Consider  the  set  of  all  possible  permutations  of  the  elements  of  all  possible 
lb-dimensional  subsets  of  2s' 

Pk  =  •  •  •  **)  :  *ii  •••»**€  Jw,  Jm  6  K*}  (4.3.1) 

Note  that  the  cardinality  of  Pk  is  equal  to  ^  ^  j  •  (k  -  1)1.  We  shall  assume  that 
the  elements  of  Pk  have  been  ordered  in  some  fashion,  and  shall  denote  these 
elements  by  n*y,  *  =  1,2,  ■  •  ■ ,  ^  ^  j ,  j  =  1,2,  •  •  • ,  (k  —  1)1.  Given  one  of  the 
subsets  Jw  €  Ms,  and  a  permutation  IIm/  of  the  elements  of  that  subset,  define 
the  corresponding  interaction  parameter  of  orderk : 

Mn*i)  =  •  •  •  aihil)1,k  (4.3.2) 


and  denote  the  largest  of  these  by 

=  max  /i(II*y)  (4.3.3) 

nwye/Tb 

Define  also  the  corresponding  loop  of  length  k: 

£(!**/)  =  •  •  •  Mikii  Ay,  (4.3.4) 

Two  loops  L{ II«y)  and  L(IIj,*n)  are  said  to  be  nontouching  if  =  $• 


52 


Note  that  there  ie  a  one-to-one  correspondence  between  the  set  of  loops  and  the 
set  of  interaction  parameters. 


Furthermore,  the  gain  in  each  loop  is  related  to  the  sise  of  the  corresponding 
interaction  parameter. 

Proposition  4.8.1:  Consider  a  permutation  I1m/  of  the  integers  JM,  and  the 
associated  loop  L(IIm/)*  Then 

•UP  *1(11*/)]  =  7V(n*/)  (4.3.5) 


Example  4.8.2  Let  ns  illustrate  our  notations  using  the  following  example. 
Consider  N  =  3,  2*  =  {1,2,3}. 

Loops  of  length  1: 

«i  =  JU  =  {»},  nw  =  (0 

L(IIui)  =  Mu&i,  m(I1ui)  =  ou,  Hi  =  ma x<o« 

Loops  of  length  2: 

H«  =  {{1,2},{1,3},{2,3}} 

Jn  =  {1»2},  Jjj  =  {1,3},  =  {2,3} 

n,u  =  (12),  nm  =  (13),  Ilasi  =  (23) 

L(IIjii)  =  — MijAjAfjxAi 

L(IIjji)  =  -Afi|A|AfjiAi 
L(Jlul)  =  — MjjAjAfjjAj 

M(n„l)  =  y/<T\lOt\ 

M(n»l)  =  y/Q\%Ol\ 

M(n»i)  =  y/onffti 
Hi  =  ma Xij  y/OijOji,  i  ^  j 
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Loops  of  length  3: 


Ks  =  {{1,2,3}},  J»i  =  {1,2,3} 

IIsii  =  (123),  n,i,  =  (132) 

L(ngn)  =  MuA)M]}A}MsiAi 
L{Um)  =  MiiA|MsiA]MhAi 
M(n#n)  =  (^n^u^ai)1/3 
M(n#ia)  =  (tfis^siffJi)1^* 

Ms  =  rnaxijj,{aiioikOkt)1,9t  ^  J  /  * 

To  illustrate  A7|/t.  and  M\jh0  consider  Jij.  Then 

Mn  0  Afis 
M\j„  =000 

M$  %  0  Mss 

and 

A1  0  0 

A,| /»  =  '  0  0  0  :  llArlU  <  7,  II A,|U  <  7  ■ 

ll  0  0  A, 

When  all  the  uncertainties  are  scalan,  then  so  are  the  blocks  of  M,  and  the 
preceding  definition  of  a  loop  reduces  to  that  used  in  Mason’s  gain  formula  (e.g., 
Theorem  0.15.20,  [ZaD63j),  which  may  be  applied  to  obtain  a  useful  expression 
for  det[/  +  MA). 

Proposition  4.S.S  Assume  that  n*  =  1,  *  =  1,2, •  •  • , N  m  (4.2.1).  Then, 

det[I  +  M  A]  =  1  +  (^(A)  +  G*{  A)  + .  •  •  +  GW(A)  (4.3.6) 


where 


G*(A)=  £  £l(ii*,)  +  **(a) 

i=l  i=i 


(4.3.7) 


and  J2*(A)  is  a  sum  of  products  of  nontouching  loops  of  length  <  k,  with  the 


property  that  each  product  contains  exactly  k  A^’s. 


Proof:  from  Mason’s  formula,  it  follows  that  det(J  +  MA)  can  be  expanded 
into  a  sum  of  products  of  nontouching  loops.  The  result  is  obtained  by  grouping 
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these  products  into  in  appropriate  G*(A),  depending  upon  how  many  A<’s  each 
product  contains. 


It  is  clear,  by  construction,  that  each  G*(A)  represents  the  *th  order  effects 
of  uncertainty  upon  det(/  4*  MA).  Furthermore,  each  G*(A)  is  a  sum  of  terms 
involving  loops  of  length  <  *.  Since  the  maximum  gain  in  each  loop  is  determined 
by  the  corresponding  interaction  parameter  (Proposition  4.3.1),  it  follows  that  if 
the  interaction  parameters  of  order  <  fc  are  small,  then  the  fcth  order  effects  will 
be  small  also.  We  summarise  this  observation  by  appealing  to  the  limiting  case. 


Proposition  4.8.4:  Let  M  and  A  be  defined  as  in  Section  4.2,  with  n<  =  1,»  = 
1, . . . ,  N.  Suppose  that 

H i=0  ,  /si,...,*  (4.3.8) 

Then,  VA  €  AM, 

G*(A)=  0,1  =  1,...,*  (4.3.9) 


Keeping  the  interaction  parameters  of  order  <  *  small  is  not,  in  general,  a 
necessary  condition  for  insuring  that  the  *th  order  effects  of  uncertainty  are  small. 
The  reason  for  this  is  that  the  preceding  analysis  takes  only  gain  information  into 
account.  The  following  example  illustrates  this  point. 

Example  4.8.5:  Consider  the  following  interconnection  matrix  for  N  =  3, 
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«*  =  !,»  =  1,2,3: 


M  = 


Oil 
-10  1 
110 


In  this  case  Ms  =  1«  hat  G*(A)  =  0.  This  example  shows  that  the  third-order 
effects  of  uncertainty  may  equal  sero  even  though  the  interaction  parameters  of 
order  three  are  nonsero. 


As  we  shall  see,  in  the  multivariable  case  directionality  as  well  as  phase  in¬ 
formation  is  ignored.  More  importantly,  our  analysis  implicitly  that  the 

Aft/a  are  mutually  independent.  In  fact,  they  are  each  functions  of  the  plant, 
compensator,  and  weightings.  The  implications  of  this  fact  will  be  discussed  in 
Section  4.5. 

Extending  Proposition  4.3.3  to  the  multivariable  case  is  possible  using  the 
formula  for  the  determinant  of  a  block  partitioned  matrix  (eg.,  [Kai80],  p.  650). 
However,  the  complexity  of  the  ensuing  expressions  renders  analysis  problematic. 
We  choose  to  adopt  an  alternate  approach,  by  redwing  the  multivariable  problem 
to  one  with  scalar  blocks.  To  do  this  we  must  invoke  a  preliminary  result. 

Lemma  4.3.6:  Consider  the  matrices  M  and  A  defined  in  Section  4.2  evaluated 
at  a  fixed  frequency.  Suppose  there  exists  A  €  AM  such  that  (i)  det[/+MA]  =  0 
and  (ii)  9[A]  =  a.  Then  there  exists  A  6  Aw  satisfying  (i)  and  (ii)  with  the 
additional  property  that  &  =  diagA*,  where  A<  =  y<xf  6<,  where  4  €  C ,  e 
d*.  and  IlftHi  =  Uzfll,  =  1. 

Proof:  Follows  from  the  results  in  Section  4  of  [Doy82j.  ■ 
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One  consequence  of  the  above  lemma  is  that  the  definition  of  the  structured 
singular  value  in  Section  4.2  may  be  modified  by  replacing  the  minimisation  over 
A  €  A*  with  a  minimisation  over  and  yt.  Define  modified  interconnection 
and  uncertainty  matrices  by 

JJ  =  {Mij  =  jfMiiVi  :  i,j  =  1  (4.3.10) 

and 

S  =  diagfa/*,  :i=l,...,N}  (4.3.11) 

Corollary  4.8.7:  Consider  M  and  A  defined  in  Section  4.2.  Then 

min  {fr[A] :  det [I  +  M A]  =  0}  =  min  {9(5] :  det[J  +  AfA]  =  0}  (4.3.12) 


It  follows  from  (4.3.12)  that  Proposition  4.3.4  may  be  extended  to  multivari¬ 
able  systems.  Specifically,  we  may  apply  Proposition  4.3.3  to  det[J  +  A/ A],  and 
define  <5*  (5)  analogously  with  (4.3.7). 

Proposition  4.8.8:  Consider  M  and  A  defined  in  Section  4.2,  the  corresponding 
interaction  parameters  (4.3.2-3),  and  the  matrices  A7  and  S  defined  by  (4.3.10- 
11).  Suppose  that 

ft=0,  1  =  1,...,*  (4.3.13) 

Then,  V&  €  G ,  t  =  l,...,J7, 

S*(S)  =0  1  =  1,...,*  (4.3.14) 
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Proof:  Each  7*  (2)  is  the  sum  of  products  of  nontouching  loops  of  order  <  l. 
Since  |M,-|  <  P(My],  it  follows  from  Proposition  4.3.1  that  (4.3.13)  implies 
(4.3.14). 


Let  us  now  derive  a  result  which  emphasises  the  difficulty  inherent  in  ana¬ 
lyzing  robustness  for  a  system  with  a  large  number  of  uncertainties.  Consider 
an  analysis  procedure  whereby  we  analyse  robustness  against  uncertainties  taken 
one  at  a  time,  then  two  at  a  time,  and  so  forth.  At  the  kth  step  of  this  proce¬ 
dure,  one  might  find  that  the  system  is  robustly  stable  against  all  combinations 
of  /ewer  than  k  uncertainties,  but  may  be  destabilised  by  some  combinations  of 
k  small  uncertainties.  The  following  result  describes  the  extreme  case  of  this 
phenomenon. 


Proposition  4.8.9:  Given  a  permutation  Iljwy  =  (»j*j  •  •  •  *jy)  6  Pn,  consider  the 
set  of  interconnection  matrices 


M(nwy)  =  | 


Mm,  T6  OjMi,*,  J*  0,. 
Mlm  =  0,  otherwise 


/  o 


(4.3.15) 


Then,  for  all  M  €  <M(IIjwy), 


Pk  —  0,  k  —  1, 2, . . . ,  N  1 


(4.3.16) 


and 


Us  =  /*(!!*</)  >  0 


(4.3.17) 
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Furthermore, 


t 


m[M|/w]  =0,  *  =  1,...,JV-1,*  =  1,... 


(4.3.18) 


and 

M(M1  =  Vn  (4.3.19) 


It  follows  from  (4.3.16-17)  and  our  earlier  discussion  that  the  1st,  2nd, . . .  (N— 
l)st  effects  of  uncertainty  upon  det(J  +  MA)  are  equal  to  sero.  Furthermore, 
stated  in  words,  (4.3.18-19)  imply  that  a  system  whose  interconnection  matrix 
satisfies  (4.3.15)  cannot  be  destabilised  by  any  combination  of  fewer  than  all  N 
uncertainties. 

Proof:  It  follows  easily  from  (4.3.15)  and  the  definition  of  the  nu*  that  (4.3.16- 
17)  are  satisfied.  It  is  also  straightforward  to  verify  that  (4.3.18)  holds.  Rom 
Proposition  4.3.8,  we  see  that 

det[7  +  A/A]  =  1  +  5^  (A) 

Selecting  the  xfn  and  y<’ s  so  that  =  0  [  Af<  ,  A7fjr<l  =  ff[A/<w<1]  and 

setting  ^  =  (pflfyry))1/*,  »  *  1,...,AT  -  1  and  SN  =  -(/ifn^y))1/*  yields 
det[J  +  A/A]  =  0.  Hence  y*(A/]  >  yt(n^y).  But  it  is  easy  to  verify  that  if 
max*  |4 1  <  l//t(IIjwy),  then  det[/-f- A/A]  ^  0.  Hence  (4.3.19)  is  satisfied. 
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4.4  Bounds  on  the  Structured  Singular  Value 


In  thia  section  we  derive  bounds  upon  the  structured  singular  value  defined  in 
(4.2.2).  We  first  note  that  most  algorithms  (e.g.,  [FaT86],  [FaT88],  [Hel88])  for 
computing  the  structured  singular  value  actually  compute  the  upper  bound  in 
(4.2.5),  which  Doyle  [Doy82]  has  conjectured  is  within  15%  of  the  true  value  of 
p\M\.  Hence,  instead  we  derive  bounds  upon  the  quantity  infers  9[DMD~X]  and 
argue  that  the  qualitative  information  furnished  by  these  bounds  remain  useful 
in  our  methodology. 

Proposition  4.4.1  Suppose  that  the  infimum  in  (4.B.5)  is  achieved.  Then, 

mf  9[DMD~X]  >  max  {  Mi,  Ms*  ■  •  ■ » Mn  }  (4-4-l) 

Proof:  First  note  that  inf^p  0\DMD~x]  =  d[t)Mt)~x\  for  some  t)  =  diag[dilni, 

•  •  • ,  d*/„w]  €  D,  if  the  infimum  is  achieved.  It  then  follows  that  inf  DeP  9[DMD~X)  > 
(di/dj)oij,  Vi,j.  Hence 

(  D€D  9[DMD~X\  )  >  {dijdijffiii,  •  (4,M,)^a<,  •  •  *  {dijdijo^i, 

=  /**(n *,) 

VTlkij  €  Pit-  This  completes  the  proof. 


We  note  that  this  proof  technique  consists  of  finding  a  lower  bound  upon  the 
quantity  inf^e?  ff\DMD~x},  and  thus  results  in  lower  bounds  upon  the  structured 
singular  value  for  the  case  of  N  <  3  (cf.  Chapter  2).  It  follows  that  these  lower 
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bounds  can  be  used  in  analysing  the  effects  of  any  two  or  three  uncertainties  in 
the  structured  singular  value  problem  for  N  >  3.  Specifically, 


m(M]  >  max  {«ll> *•><(•<} 


(4.4.2) 


Alternately,  one  can  also  use  (4.2.6)  to  derive  lower  bounds  upon  the  structured 
singular  value,  as  done  in  [Fre89a]  for  the  case  N  —  2.  However,  the  bounds  ob¬ 
tained  in  this  way  are  sometimes  complicated  and  thus  may  obscure  information 
we  desire  from  the  interaction  parameters. 

Proposition  4.4.2  Suppose  pn  >  0.  Then 

<  „„  £ :(^)‘  (4.4.3) 

Furthermore,  if  Pk/Ps  <  1,  VJfc,  then 

P[M\<Y,Pk  (4.4.4) 

k=i 


Proof:  See  Appendix  G. 


Proposition  4.4.2  provides  an  upper  bound  useful  in  the  case  that  the  most 
significant  effects  are  due  to  the  interaction  of  all  AT  uncertainties,  i.e.,  pu  < 
Ps,  Vh  =  1, 2,  •  •  • ,  N  —  1.  Bounds  useful  in  other  cases  may  also  be  derived  by 
following  a  procedure  suggested  in  Chapter  2.  Note,  however,  that  the  bounds 
obtained  in  such  cases  may  involve  significant  compl  exity,  especially  when  the 
number  of  uncertainties  is  large.  To  summarise,  it  follows  from  these  lower  and 
upper  bounds  that  large  values  of  the  interaction  parameters  can  lead  to  a  large 
structured  singular  value,  and  thus  robustness  problems.  This,  indeed,  is  consis¬ 
tent  with  our  analysis  in  the  preceding  section.  Consequently,  the  bounds  derived 
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in  this  section  along  with  the  analysis  in  the  preceding  section  do  show  that  the 
interaction  parameters  can  qualitatively  characterize  robustness  properties.  A 
useful  application  is  to  analyze  these  parameters  to  gain  design  insights  into  a 
given  problem,  as  in  [Fre89a] ,  [ChF89a],  and  [ChF89b],  In  applying  this  method¬ 
ology,  it  is  significant  to  note  that  the  blocks  M,y  of  the  interconnection  matrix 
will  be  mutually  interrelated,  and,  as  in  the  cited  references,  it  is  important  to 
analyze  this  interrelation  to  obtain  design  insights. 

4.5  An  Analysis  Example 

We  now  illustrate  our  analysis  procedure  fay  studying  a  three-block  structured 
singular  value  problem.  Consider  the  linear  time-invariant  feedback  system  de¬ 
picted  in  Figure  3.1.  The  transfer  functions  P(s)  and  F(s)  are  those  of  the  plant 
model  and  compensator,  respectively.  We  shall  assume8  that  the  plant  is  square, 
and  that  detP(s)  £  0.  The  signals  r(«),  d(a),  and  y(a)  are  the  reference  input, 
disturbance  input,  and  system  output,  respectively.  Define  the  input  open  loop 
transfer  function,  sensitivity  function,  and  complementary  sensitivity  function: 

£,(.)  =  F(«)P(»),  S,M  =  |/  +  £,(»)]-',  IY(.)  =  £,(.)]/  +  £,(«)]-'  (4.5.1) 

and  the  output  open  loop  transfer  function,  sensitivity  function,  and  complemen¬ 
tary  sensitivity  function: 

Lo{b)  =  P(s)F(s),  So{»)  =  [I  +  Lois))-1,  T0{s)  =  L0{*)[I  +  L0{s)]~l  (4.5.2) 

sOur  results  extend  to  nonsquare  plants  that  satisfy  the  relevant  assumption  of  left  or  right 
invertibility. 
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We  consider  a  benchmark  robust  performance  problem  studied  in  [Fke89a], 
[Fre89bj,  [SMD88]  and  [Ste87]  with  an  additional  source  of  modelling  error  which 
we  take  to  be  an  unstructured  additive  uncertainty.  As  discussed  in  [DWS82], 
such  a  performance  robustness  problem  can  be  treated  as  an  equivalent  robust 
stability  problem  against  the  actual  sources  of  plant  uncertainty  plus  an  addi¬ 
tional  source  of  fictitious  uncertainty  used  to  represent  the  performance  spec¬ 
ification.  In  thus  regard,  the  present  problem  is  equivalent  to  the  problem  of 
maintaining  stability  robust  against  plant  uncertainty  in  the  form  of 

Pt(s)  =  [/  +  r,(s)  Aj(s)]-1[P(s)  +  r,(s)  A*(«)][/  +  n(a)  Ax(«)]  (4.5.3) 

which  corresponds  to  the  block  diagram  depicted  in  Figure  4.1.  Here  r<(«), 
i  =  1,2,3,  are  stable,  minimum  phase  scalar  weighting  functions  used  to  de¬ 
scribe  how  the  uncertainties  and  performance  requirement  vary  with  frequency, 
and  A<(s),  t  =  1,2,3  are  unstructured  uncertainties  which  satisfy  ||A*||oo  <  7. 
Note  that  Ai(«)  and  As(«)  are  actual  sources  of  plant  uncertainty,  while  At(s) 
is  the  fictitious  uncertainty  used  to  represent  the  performance  specification.  Our 
nominal  performance  goal  will  be  achieved  if  the  output  sensitivity  function  sat¬ 
isfies  a  frequency  dependent  bound 


IMolioo  <  1  h  (4-5-4) 

and  the  robust  performance  goal  will  be  satisfied  if  VAi,  As,  || Ai||oo  <  7* 
ilAsIloo  <  7» 

IksSolloo  <  1/7  (4.5.5) 

with  S'0  =  [I  +  FF]-1  and  F  =  [P  +  r*A»l[ J  +  r^j] 
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We  shall  now  discuss  the  stability  robustness  problem  in  detail.  For  purposes 
of  illustration,  we  shall  mainly  demonstrate  this  procedure  which  centers  at  an¬ 
alysing  the  interaction  parameters;  a  complete  derivation  of  design  insights  for 
this  robust  performance  problem  is  found  in  [ChF89b]. 


We  now  formulate  this  stability  robustness  problem  as  a  three-block  struc¬ 
tured  singular  value  problem.  Our  first  task  is  to  find  the  corresponding  inter¬ 
connection  matrix.  Manipulating  the  block  diagram  in  Figure  4.1  yields 


M  = 


t\Tj  -riTjP-1  riT/P-1  ' 
—r%SoP  rjS0  -fjSo 
-r,5/  -rsTiP-1  raTrP~l . 


(4.5.6) 


Next,  we  need  to  identify  the  interaction  parameters  defined  in  (4.3.2-3),  and  for 
convenience  we  shall  adopt  the  ordering  of  these  parameters  given  as  in  Example 
4.3.2: 

Mi  =  max  {<T[ri7>],  *[r*So],  ^rsT/P-1]  > 

<*(n,u)  =  y'a[r,S<,P|i>[r1r,P-i] 

m(IIui)  =  V»[f.g,lJ|r1r,i’-1|  .  . 

M(n»,)  =  v'«[r,Sb]»(r,r,P-,l  1  ■  '  ’ 

M(n.u)  =  (»|i-ir,p-1|»(r,s0|»[r,s,|)l/* 

n(ns„)  =  (»|r1r,p-‘|ff(r,s0p]«>[r1r,p->)),/* 


We  first  analyse  the  effects  that  are  due  to  each  individual  uncertainty.  Note 
that  if  the  feedback  system  is  nominally  stable,  then  it  will  be  robustly  stable 
against  Aj  alone  if  and  only  if  ||riT/||oo  <  1/7  [DoS81].  Similar  statements  apply 
to  A]  and  ||rjS0||oo>  and  to  As  and  ||rsT/.P-1||00.  Hence  summarises  stability 
robustness  with  respect  to  each  of  the  three  uncertainty  sources  individually. 


Next,  we  analyse  the  effects  of  uncertainties  taken  two  at  a  time.  The  effect 
of  interactions  between  Ai  and  As  is  quantified  by  n[M\jtl]  and  has  been  studied 
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in  [Fre89a],  [Fre89b].  Applying  the  bound  (4.4.2-3)  for  N  =  2  to  this  problem,  it 
follows  that  the  parameter  /i(IIsii)  determines  the  extent  to  which  these  two  un¬ 
certainties  can  interact  to  cause  stability  robustness  difficulties  even  if  the  system 
is  robustly  stable  against  Ax  and  A*  individually.  To  prevent  this  undesirable 
effect  imposes  a  nontrivial  problem,  for  merely  requiring  9[r\Ti]  and  o[rtSo\  to 
be  small  does  not  imply  that  fi(ll)u)  is  also  small.  To  this  point,  let  the  plant 
condition  number  [GoV83]  be  defined  by  ic[P(;(t>)]  =  9[P{juf)\f o{P{ju)\. 
Proposition  4.5.1  [Fre89a]  Let  /i(nlu)  be  defined  in  (4.5.7).  Then,  at  each 
frequency, 

M(n,u)  <  yfiW  •  (4.5.8) 

Since  the  bound  (4.5.8)  may  be  tight  [Fre89a],  it  follows  that  stability  robustness 
against  Ax  and  As  together  is  potentially  much  worse  than  that  against  each 
of  these  two  individually  at  frequencies  for  which  the  plant  is  ill-conditioned 
(w[P]  >  1).  It  was  shown  in  [Fre89a],  [Fre89b]  that  requiring  p(IIsXX)  to  be  small 
imposes  more  stringent  constraints  upon  the  closed-loop  transfer  functions  than 
those  that  are  due  to  the  stability  robustness  requirements  with  respect  to  Ax  or 
Aa  alone  (i.e.,  ||rxr/||oo  <  1/7,  ||rfe3o||a»  <  1/7.  respectively). 

Similarly,  we  analyse  the  combined  effects  of  uncertainties  Aj  and  As. 
Proposition  4.5.2  Let  p( IIssi).  Mi  defined  in  (4.5. 7).  Then,  at  each  fre¬ 
quency, 

M[M\jJ  <  2pi  (4.5.9) 
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Proof:  By  definition, 


/*( n»i)  =  \Jo\riSo\o\rtTjP-1]  < 
The  proof  now  follows  from  (4.4.3)  for  N  =  2. 


Hence,  Proposition  4.5.2  reveals  that,  »n1iW  the  problem  of  maintaining  sta¬ 
bility  robustness  with  respect  to  At  and  At,  a  feedback  design  which  is  robust 
against  As  and  A*  alone  will  tend  to  be  robust  with  respect  to  these  two  uncer¬ 
tainties  taken  together,  independent  of  the  sise  of  the  plant  condition  number. 

Lastly,  we  consider  the  combined  effects  of  uncertainties  Ai  and  A*. 
Proposition  4.5.3  Let  |i(nni),  be  defined  in  (4.5. 7).  Then,  at  eaeh  fre¬ 
quency, 

<  ^»i(l  ri  I  +<*i)  (44.10) 

x(nMI)  >  (I  r,  I  -B[r,r,])  (44.11) 


Proof:  It  follows  from  definition  that 

M(nm)  =  \J 9{riTjP-l\o[rtSi] 

The  proof  is  now  completed  by  applying  the  triangle  inequality. 


The  bounds  in  (4.5.10-11)  exhibit  that  stability  robustness  with  respect  to 
At  and  Aa  together  is  potentially  much  worse  than  that  with  respect  to  either  of 
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these  uncertainties  alone  at  any  frequency  for  which  |  ri(/w)  |»  1/7.  As  shown 
in  [ChF89b],  maintaining  /t(nju)  small  at  such  a  frequency  imposes  a  more 
stringent  constraint  than  what  is  necessary  for  achieving  stability  robustness 
against  As  alone  (UraTjjP- 1||00  <  1/7).  However,  it  was  shown  in  [ChF89b  that 
the  two  parameters  m(Hisi)  and  /i(IIiii)  are  interrelated  in  such  a  way  that 
keeping  both  of  them  small  at  such  a  frequency  imposes  similar  requirements 
upon  the  closed-loop  transfer  functions,  and  thus  maintaining  A*(n«i)  small  is 
also  useful  for  keeping  /x(IImi)  small.  Hence,  achieving  stability  robustness  with 
respect  to  Aj  and  As  together  also  tends  to  guarantee  a  certain  smaller  level  of 
stability  robustness  with  respect  to  the  combined  effects  of  Ai  and  As. 

Finally,  we  need  to  consider  the  simultaneous  effects  of  all  three  uncertainties. 
Recall  that  these  effects  are  determined  by  the  parameters  p(IIsii)  and  p(Hsu)‘ 
Fortunately,  for  this  particular  problem,  requiring  Hi  and  £4  to  be  small  forces 
Hi  to  be  small  as  well. 

Proposition  4.5.4  Let  Hu  H*>  and  Ps  te  defined  in  (4.8.8).  Then,  at  each 
frequency, 

H»  <  max  {hi,  Hi  }  (4.5.12) 


Proof:  First,  it  is  straightforward  to  show 

M(n$«)  <max  {HuPi} 


since 


M(n,„)  =  (elr1r,p-'|»IrBT,p-1|«[r1So/,])1/' 
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Next,  following  the  same  line  establishes 

M(n#n)  <  max  {/ii,Ms> 

since 

=  (^(nuiJ^lrjSo])  ^ 

<  (m5  -Mi)1/S 


We  now  close  this  section  by  summing  up  the  analysis  procedure  demon¬ 
strated  as  above.  Essentially,  we  decompose  a  complicated  robust  performance 
problem  into  a  sequence  of  more  tractable  problems  in  which  we  analyse  the 
interaction  parameters.  An  important  step  of  this  method  is  to  study  interrela¬ 
tions  that  exist  among  all  interaction  parameters  and  thus  to  obtain  information 
useful  for  purposes  of  design.  For  a  particular  robustness  problem,  interrelations 
between  blocks  of  the  interconnection  matrix  do  allow  us  to  draw  significant  in¬ 
sights  into  the  problem.  Finally,  we  note  that  the  analysis  in  this  section  also 
echoes  Example  4.3.5  in  showing  why  keeping  interaction  parameters  of  order  <  k 
small  is  not  in  general  necessary  for  keeping  the  fcth  order  effects  of  uncertainties 
small. 


60 


4.6  Concluding  Remarks 


We  have  developed  an  analysis  method  for  analyzing  robustness  problems  asso¬ 
ciated  with  multiple  sources  of  modelling  error.  We  used  Mason's  gun  rule  to 
identify  a  set  of  interaction  parameters  which  were  shown  to  capture  the  essen¬ 
tials  of  such  problems.  These  parameters  were  further  used  to  derive  bounds 
upon  the  structured  singular  value  and  their  importance  was  manifested  through 
analyses  of  certain  extreme  cases.  The  analysis  methodology  we  propose  is  to 
decompose  a  robustness  problem  into  a  sequence  of  problems  in  which  we  ana¬ 
lyze  the  interaction  parameters.  An  important  step  of  this  method  is  to  study 
interrelations  among  all  interaction  parameters  and  thus  draw  insights  useful  in 
design.  A  benchmark  robust  performance  problem  was  analysed  to  illustrate  this 
method. 
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Chapter  5 

Achieving  Robust  Performance 
via  H°°/H2  Mixed  Sensitivity 
Optimization 

5.1  Introduction 

In  this  chapter  we  shall  consider  the  robust  performance  problem  of  achiev¬ 
ing  good  command  following  (equivalently,  small  output  sensitivity)  despite  un¬ 
structured  multiplicative  uncertainty  at  the  plant  input.  This  is  a  difficult  design 
problem,  because  achieving  robust  performance  requires  that  feedback  properties 
at  both  plant  input  and  output  be  manipulated  simultaneously.  This  precludes 
the  naive  use  of  multivariable  loop-shaping  techniques  as  well  as  of  synthesis 
techniques  such  as  H*/H°°  mixed  sensitivity  optimisation.  The  goal  of  robust 
performance  may  be  translated  into  one  of  minimising  the  structured  singular 
value  (SSV)  [Doy85].  However,  direct  minimisation  of  the  SSV  is  still  an  exper¬ 
imental  procedure,  with  many  difficulties  remaining  to  be  resolved  [DLP86]. 

In  [Re89a],  [Fre89b],  and  [Fre90],  a  multivariable  loop-ehaping  approach  to 


71 


thin  problem  was  presented,  based  upon  an  earlier  version  of  the  bounds  derived 
in  Chapter  2.  The  approach  taken  was  to  analyse  the  interrelations  among 
the  blocks  of  the  interconnection  matrix  to  obtain  guidelines  for  compensation 
selection.  In  this  chapter  we  shall  adopt  an  alternate  approach,  and  analyse 
the  bounds  to  obtain  guidelines  for  selecting  weighting  functions  to  be  used  in 
H*/H°°  mixed  sensitivity  optimisation.  These  weighting  functions  will  depend 
upon  the  original  weighting  functions  used  to  represent  design  specifications  in 
the  robust  performance  problem,  as  well  as  the  directional  properties  of  the  plant. 
We  shall  show  that  tnitiimising  the  H*  or  H°°  norm  of  the  mixed  sensitivity 
function  with  these  new  weightings  will  lend  to  minimise  the  sise  of  the  structured 
singular  value  for  the  original  robust  performance  problem.  Results  summarised 
in  this  chapter  appear  in  [L0F88],  [LoF89a],  (LoF89b}> 

5.2  Simultaneous  Uncertainty  and  Robust  Per¬ 
formance  Problems 

In  this  section  we  review  relevant  results  from  [Fre89a]  and  [Fre89b].  The  math¬ 
ematical  problem  we  study  has  two  distinct  physical  interpretations,  either  as 
a  simultaneous  uncertainty  problem  or  as  a  robust  performance  problem.  In 
the  former,  we  ask  that  the  Bystem  be  robustly  stable  against  modelling  errors 
occurring  at  two  points  in  the  feedback  loop.  In  the  latter  we  ask  that  a  perfor¬ 
mance  property  defined  at  one  point  be  robust  against  uncertainty  occurring  at 
the  other.  Both  interpretations  are  useful  in  guiding  and  explaining  subsequent 
developments. 
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Our  first  task  is  to  define  several  important  transfer  functions.  Consider 
the  linear  time-invariant  feedback  system  pictured  in  Figure  5.1.  The  transfer 
functions  P(s)  and  K(s)  are  those  of  the  plant  model  and  compensator, 
respectively.  We  shall  assume1  that  the  plant  has  n  inputs  and  outputs,  and  that 
det  P(s)  0.  The  signals  r(s),  e(s),  and  y(s)  are  the  reference  input,  error 
signal,  and  system  output,  respectively.  Properties  of  this  system  are  governed 
by  two  sets  of  transfer  functions.  Breaking  the  loop  at  the  plant  input  yields 
the  input  open-loop  transfer  function,  sensitivity  function,  and  complementary 
sensitivity  function 


£,(.)  =  K(,)P(,).  S,(.)  &  [1  +  £,(*))-*,  !>(«)  =  £,(.)  [/  +  L,(.)r‘  (6.2.1) 

Breaking  the  loop  at  the  plant  output  yields  the  output  open  loop  transfer  func¬ 
tion,  sensitivity  function,  and  complementary  sensitivity  function 

LoW  =  P(«)  *(.),  So(.)  =  (/+£oW)-‘,  r„(.)  =  £c(.)|/+L<,(.))-1  (6.2.2) 

We  shall  suppose  that  the  true  plant  differs  from  the  model  due  to  uncertainty 
that  is  present  simultaneously  at  the  plant  input  and  plant  output.  Specifically, 
we  suppose  that 

P*(«)  =  (/  +  r,(,)A,W]-‘  P(.)[/  +  r,  W  A,(.)]  (6.2.3) 

A  block  diagram  of  the  associated  feedback  system  is  in  Figure  5.2.  The  transfer 

functions  Ai(s)  and  Aj(«)  used  to  model  system  uncertainty  will  be  referred  to 

1Oar  results  extend  to  nonsqnsre  plants  that  satisfy  the  relevant  assumption  of  left  or  right 
invert  ibilify. 
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Figure  5.1:  Block  Diagram  of  Feedback  System 
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Figure  5.2:  Block  Diagram  With  Uncertainty 
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as  input  and  output  perturbations,  respectively,  and  are  assumed  to  be  stable 
and  proper.  In  general,  the  size  of  each  uncertainty  will  vary  with  frequency. 
This  property  is  modelled  by  using  the  weighting  functions  r<(«).  With  no  low* 
of  generality,  these  functions  are  assumed  to  be  stable  and  minimum  phase, 
although  they  may  be  improper  [DWS82].  The  sources  of  each  uncertainty  are 
now  described. 

Uncertainty  of  the  form 

Pi(«)  =  P(a)  [I  +  ri(s) Ai(e)]  (5.2.4) 

is  termed  multiplicative  input  uncertainty.  It  typically  arises  because  of  unmod¬ 
elled  actuator  dynamics  or  high  frequency  plant  modelling  errors.  Suppose  that 
the  feedback  system  is  nominally  stable  and  that  the  perturbation  Aj  may  be  any 
stable  transfer  function  satisfying  the  upper  bound  fr[Ai(ju;)]  <  7,  Vw.  Then  the 
system  will  be  robustly  stable  if  and  only  if  the  input  complementary  sensitivity 
function  satisfies  the  bound  [DWS82]: 

9[ri{jv)TiO'v)]  <  l/l  ,  Vw  (5.2.5) 

Hence  6[fi(jw)^V(iw)]  is  inversely  proportional  to  a  frequency  dependent  stabil¬ 
ity  margin  against  uncertainty  (5.2.4) 

Uncertainty  of  the  form 

P,  (*)  =  [/  +  r,M  A,(.)j-‘/>(.)  (5.2.6) 

is  termed  divisive  output  uncertainty.  Its  presence  may  be  due  to  modelling 
errors  at  the  plant  output.  For  our  purposes,  however,  it  is  used  to  represent  a 
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performance  specification  such  as  an  upper  bound  upon  the  response  of  the  error 
signal  to  the  reference  signal  shown  in  Figure  5.1  [DWS82].  Suppose  that  the 
feedback  system  is  nominally  stable  and  that  the  perturbation  Ai  may  be  any 
stable  transfer  function  satisfying  the  upper  bound  &[Aj(j'w)]  <  7,  Vw.  Then 
the  system  will  be  robustly  stable  if  and  only  if  the  output  sensitivity  function 
satisfies  the  bound  [DWS82]: 

0[rt{ju)S0[ju>)  <  I/7  ,  Vw  (5.2.7) 

Hence  8[ri{ju)So(Ju)\  is  inversely  proportional  to  a  frequency  dependent  sta¬ 
bility  margin  against  uncertainty  (5.2.6).  Alternately,  if  the  output  uncertainty 
arises  from  an  attempt  to  represent  a  performance  goal,  we  say  that  the  goal  of 
nominal  command  following  is  satisfied  if  (5.2.7)  holds. 

Testing  for  robust  performance  requires  use  of  the  structured  singular  value 
[DWS82].  Rearranging  the  block  diagram  in  Figure  5.2  into  the  form  shown  in 
Figure  2.1,  we  have1 


rxTj  —t\TiP~1 

—t%SqP  f|  Sq 


(5.2.8) 


and 

A  =diag[A,,A,]  (5.2.9) 

Defining  the  structured  singular  value  as  in  Section  2.2,  we  have  that  fr[r,50]  < 
1/7,  VAi  such  that  9(Ai]  <  7  if  and  only  if  p[M]  <  1/7.  Hence  the  dmign 
goal  of  achieving  robust  performance  can  be  translated  into  one  of  minimising 
3  Am  ami,  we  suppress  dependence  on  frequency  where  convenient. 
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Moo  =  supw  n[M(jw)}.  We  shall  approach  this  problem  by  analysing  the  bounds 
upon  the  structured  singular  value  presented  in  Chapter  2. 


Lemma  5.2.1:  With  M  and  A  defined  by  (5. 2. 8-9),  the  structured  singular 
value  satisfies  the  bounds 


and 


* 

V 

n[M]  >  max  • 


9[rI2>|,9|r,£,]  , 
v/9[r1r,J>-*]9[r,S,l>] 


1. 


(5.2.10) 


/*[Af]  <  y/slnTjP-^OlrjSoP] 
+  max{ff  [riT/],  [rj5o]} 


(5.2.11) 


It  follows  from  Corollary  5.2.1  that  if 

ffl^TjP'^lriSoP]  >  max{ff[rirj],^[rj5o]} 


and 

OfrMP-'WriSoP]  >  1  (5.2.12) 

then  performance  robustness  will  be  poor  even  if  nominal  performance  and  robust 
stability  are  satisfactory.  Hence  it  is  desirable  to  know  how  a  compensator  should 
be  selected  to  prevent  (5.2.12)  from  being  satisfied.  We  approach  this  problem 
by  deriving  a  set  of  upper  and  lower  bounds  that  are  analogous  to  (5.2.10-11)  and 
that  display  the  plant  directionality  properties.  For  simplicity  of  exposition,  we 
shall  restrict  our  discussion  to  systems  with  two  inputs  and  outputs.  Extensions 
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to  n  inputs  and  outputs  are  available.  The  following  results  are  similar  in  spirit 
to  those  of  Section  2,  [Pre89b], 

Denote  the  plant  singular  value  decomposition  by 

P  =  WTZa 

=  Eteisfr*  (5-2-13) 

i=l 

where  T  =  diag[ri,rj]  contains  the  singular  values  with  usual  ordering,  the 
columns  of  W  =  [ti^toj]  are  the  left  singular  vectors,  and  those  of  Z  =  [zi*s] 
are  the  right  singular  vectors.  Define  the  coupling  coefficients 

B(M  =  MSowtll  •  tr,«n  (SAM) 


The  following  result  follows  similarly  to  Theorem  2.2  of  [Fre89b]. 


Theorem  5.2.2:  For  each  pur  of  plant  singular  values,  the  structured  singular 
value  satisfies 


^[»’ir/],^[rjSo] 


n[M]  >  max 


(5.2.15) 


Suppose  that  we  desire  the  structured  singular  value  to  satisfy  <  1.  Then 
each  lower  bound  implied  by  (5.2.15)  corresponds  to  a  necessary  condition  for 
robustness.  By  deriving  upper  bounds  upon  the  structured  singular  value,  it  is 
also  possible  to  state  sufficient  conditions  in  terms  of  the  coupling  coefficients. 

Theorem  5.2.8:  For  each  pair  of  plant  singular  values,  the  structured  singular 
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value  must  satisfy 


**[«]< 


N 


££a(M)fa/n)  + 


k=l  J=1 


max 


(5.2.16) 


A  tighter  upper  bound  may  be  obtained  if  an  additional  orthogonality  con¬ 
dition  ia  satisfied. 


Corollary  5.2.4:  Suppose  that 


(i)  S0v)k±S0wi  ,  k  ±  l 

(ii)  TiZk±TiZi  ,  k^l 


(5.2.17) 


Then 


H[M\  <  mt*y/B{k,l){Th/Ti)  +  ““{  »MW  }  (018) 


Together,  Theorems  S.2.2-3  and  Corollary  5.2.4  show  that  requiring  cou¬ 
pling  coefficient  to  be  suitably  small  tends  to  insure  that  performance  robustness 
will  be  satisfactory  provided  that  the  bounds  ff[ri7>]  <  1  and  fr[rjSo]  <  1  are 
satisfied. 

Note  first  that  each  coupling  coefficient  is  the  product  of  the  weightings  and 
two  terms  that  are  each  the  gun  of  a  closed-loop  transfer  function  in  a  particular 
direction.  It  follows  that  the  weighted  condition  number 

#Cr[P]  =  |rir,|.*[P]  (5.2.10) 

is  a  particularly  relevant  parameter  of  the  design  problem.  At  frequencies  for 
which  (5.2.19)  is  large,  it  is  necessary  that  the  sensitivity  and  complementary 
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sensitivity  functions  be  shaped  carefully  to  prevent  the  structured  singular  value 
from  becoming  large.  This  imposes  constraints  upon  the  open-loop  singular 
values. 

Lemma  5.2.5:  For  each  pur  of  singular  values,  the  coupling  coefficients  are 
bounded  above  and  below  by 

|rir,|Sf[Lo]/|l  >  B(fc,l)  >  |rlr,|ff[L0]/(l  +  ^\L0])i  (5.2.20) 

|rjrj|?r(2//]/|l  —  2[I»/]|J  >  B (fc,l)  >  jriri|ff(L/j/(l  +  *[£/])*  (5.2.21) 

It  also  follows  from  the  preceding  theorems  that  the  closed  loop  transfer  func¬ 
tions  must  satisfy  certain  directionality  properties.  One  approach  to  achieving 
these  properties  is  to  translate  the  closed-loop  properties  into  specifications  upon 
the  open  loop  transfer  function.  This  is  explored  in  [Fre89b]  and  [FreOO].  In  the 
following  section  we  shall  pursue  an  alternate  approach. 

5.3  Robust  Performance  via  Mixed  Sensitivity 
Optimization 

We  saw  in  the  preceding  section  that  the  goal  of  robust  performance  may  be 
translated  into  one  of  manipulating  the  sise  of  the  coupling  coefficients.  In  this 
section  we  shall  show  how  this  may  be  done  using  H*/H°°  mixed  sensitivity 
optimisation  techniques.  For  simplicity  of  exposition,  we  shall  consider  only 
systems  with  two  inputs  and  two  outputs.  In  this  case,  Theorem  5.2.2  reduces 
as  follows. 


81 


Corollary  5.S.1:  Necessary3  conditions  for  achieving 


p[M]  <  1  (5.3.1) 

are  that 

max{ff[rilj],  ff[fj5o]}  <  1  (5.3.2) 

lkiT/*j||  •  ||r150tt,1|j(/c)  <  1  (5.3.3) 

Ikil/zill  •  ||r#50to,||(l/ic)  <  1  (5-3-4) 

||riT/Zij|  ■  ||rs5o«V|||  <  1  (5.3.5) 

WriTjZtW  •  || rrfoiOiW  <  1  (5.3.6) 


Note  first  that  if  (5.3.2)  is  satisfied,  then  the  only  one  of  conditions  (5.3.3-6) 
remaining  of  interest  is  (5.3.3) 

We  now  briefly  review  the  mixed  sensitivity  minimisation  problem.  Consider 
the  system  depicted  in  Figure  5.3.  The  transfer  function  from  the  disturbance  d 
to  the  weighted  control  signal  yi,  and  system  output  y>,  is  given  by 


Vi(«) 

»w 


=  G{s)d{s) 


where 


(5.3.7) 


9R«c*n  tkat  Thtonm  5.3.2  and  Corollary  5.5.4  imply  (hat  makiaf  (5.3.24)  amall  also  taada 
to  iaaara  (hat  (5.3.1)  ia  aatiafiad. 
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The  weighting  functions  Wx,Wt,  and  W4  are  assumed  to  be  stable  and  minimum 
phase.  We  shall  refer  to  (5.3.7)  as  the  weighted  mixed  sensitivity  function.  Define 
the  H°°  mixed  sensitivity  problem  as  that  of  minimising 

||G||«  =8upfr[G(;w)]  (5.3.8) 

over  all  internally  stabilising  controllers.  Similarly,  we  may  define  the  H*  mixed 
sensitivity  problem  as  that  of  minimising 

l|G||s  =  jT  trace  G^-iwJGO'w)]  1  (5.3.9) 

over  all  internally  stabilising  controllers.  Computationally  tractable  solutions 
to  these  problems  using  state-space  techniques  have  recently  become  available 
[DGK89]. 

We  shall  now  show  how  to  pick  weighting  functions  so  that  minimising  (5.3.8) 
tends  to  minise  Modifications  of  this  procedure  to  the  problem  of  minimising 
(5.3.9)  are  found  in  [LoF89b]. 

Theorem  5.8.3:  Consider  Af  («)  and  G(s)  defined  by  (5.2.8)  and  (5.3.7),  respec¬ 
tively,  and  evaluated  at  a  fixed  frequency.  Assume  that  the  weightings  in  (5.3.7) 
are  given  by 


3 

II 

* 

O  0 

►-  0 

(5.3.10) 

Wx  =  rXTt{K/a)I 

(5.3.11) 

Wt  =  r,/ 

(5.3.12) 
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and 


where  W  is  the  matrix  of  pleat  left  singular  vectors  (5.2.13),  fi  and  r*  are  the 
weighting  functions  in  (5.2.4)  end  (6.2.6),  rf  is  the  smeller  plant  singular  value, 
k  is  the  plant  condition  number,  and  *  >  |o|  >  1.  Suppose  that 

*[G]  <  1  (5.3.13) 

Then  (S.3.2-6)  are  satisfied. 


Proof:  Define  Gx  =  -WtP~lT0W4  and  Gt  =  WtS0W4.  It  is  well-known  that 


9[C]  >  maxd[G«] 

(5.3.14) 

and,  for  any  unit  vector  z, 

*(G|  >  ||C<*II 

(5.3.15) 

Notice  that 

a 

=  °/a 

(5.3.16) 

Since  ic/|a|  >  1  by  construction,  and  since  Z  is  unitary,  we  have  from  (5.3.14-15) 

*(G]  >  9[riTt] 

(5.3.17) 

and,  with  ej  =  [0  l]r, 

6(G)  >  ||Giei|| 

=  forf*(«/«)| 

(5.3.18) 

Similarly, 

$ 

II 

* 

o  Q 

>-  O 

(5.3.19) 
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which  implies  that 


9[G]  >  *(rtS0] 


(5.3.20) 


and 

0[G\  >  IMoWill  (5.3.21) 

From  (5.3.17)  and  (5.3.20),  it  follows  that  (5.3.13)  implies  (5.3.2)  and  thus  also 
(5.3.4-0).  Finally,  (5.3.13),  (5.3.18),  and  (5.3.21)  taken  together  imply  (5.3.3). 


Although  these  are  not  avffietent  conditions,  the  upper  bounds  in  Theorem 
5.2.3  suggest  that  minimMing  ||G||oo  also  tends  to  make  pm  itself  small,  and  this 
conjecture  appears  to  be  borne  out  by  our  design  examples. 

We  now  discuss  how  weighting  functions  with  the  properties  (5.3.10-12)  may 
be  constructed.  First,  it  may  not  be  possible  to  find  a  rational  matrix  Wt(a)  so 
that  (5.3.10)  holds,  because  the  singular  vectors  need  not  be  complex  analytic 
and  thus  need  not  have  frequency  responses  equal  to  those  of  rational  transfer 
functions.  However,  these  vectors  can  often  be  approximated  sufficiently  closely 
for  our  purposes  by  rational  functions.  (On  the  other  hand,  tie  smallest  singu¬ 
lar  value  rs  and  the  condition  number  k  can  always  be  approximated  arbitrarily 
closely  by  the  gain  of  a  stable,  rational,  minimum  phase  transfer  function.)  Fi¬ 
nally,  we  need  to  select  a  rational  transfer  function  a(s).  To  do  this  we  suppose 
that,  as  is  typical,  the  weighting  functions  r\  and  ri  have  the  general  shape  shown 
in  Figure  5.4.  Let  us  now  investigate  the  implications  of  requiring  (5.3.2)  to  hold 
at  a  frequency  for  which  |ri(jci/)|  >  1  or  |rj(jw)j  >  1. 
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Lemma  5.S.3:  (a)  Suppose  that  0[rs5b]  <  1  end  that  |r*|  >  1.  Then 


\\riTM  m  |n|  (5.3.22) 

(b)  Suppose  that  o[t\Tj\  <  1  and  that  |ri|  >  1.  Then 

||n*>«l||  «  M  (5.3.23) 


It  follows  that  if  we  wish  to  satisfy  (5.3.2)  and  (5.3.3)  at  low  frequencies,  for 
which  |r>|  >  1,  then  we  should  require  that  ||rs5ottf|ic||  <  1.  Hence  we  should 
have  |a|  »  k  at  low  frequencies.  Similar  considerations  show  that  we  should 
require  |a|  »  1  at  high  frequencies.  Hence  we  shall  select  a  function  a(s)  whose 
gain  decreases  from  Ja]  =  k  to  jaj  =  1. 

Of  course,  the  weightings  described  above  may  need  iteration  before  a  satis¬ 
factory  design  is  obtained. 


Finally,  we  state  dual  results  that  shall  prove  useful  in  our  example.  To  do 
this,  define  the  mixed  sensitivity  function 


G(s)  = 


-Wi(«)T,(aWa) 

Wt(a)P(s)Sj(s)W4(s) 


(5.3.24) 


corresponding  to  the  block  diagram  in  Figure  5.5. 

Theorem  5.3.4:  Consider  Af(a)  and  &($)  defined  by  (5.2.8)  and  (5.3.24),  re¬ 
spectively,  and  evaluated  at  a  fixed  frequency.  Assume  that  the  weightings  in 
(5.3.24)  are  given  by 

=  *[J  2]  («"*» 
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Figure  5.5:  Block  Diagram  for  Alternate 


WX  =  T\I 

Wx  =  (r,//Jfi)/ 


(5.3.26) 

(5.3.27) 


where  Z  is  the  matrix  of  plant  right  angular  vectors  (5.2.13),  ri  and  fs  are  the 
weighting  functions  in  (5.2.4)  and  (5.2.6),  rj  is  the  smaller  plant  singular  value, 
k  is  the  plant  condition  number,  and  k  >  \fi\  >  1.  Suppose  that 

0[&]  <  1  (5.3.28) 

Then  (S.3.2-6)  are  satisfied.  ■ 

Proof:  Similar  to  that  of  Theorem  5.3.2.  ■ 

The  transfer  function  /?(«)  should  increase  from  1  at  low  frequencies  to  ic  at 
high  frequencies. 

5.4  Design  Example  Using  H°°  Mixed  Sensitiv¬ 
ity 

We  now  apply  the  procedure  developed  in  the  preceding  section  to  the  problem 
of  designing  a  robust  performance  longitudinal  control  system  for  a  highly  ma¬ 
neuverable  aircraft.  This  example  was  originally  proposed  by  Stein  [Ste84].  The 
aircraft  characteristics  and  linearised  models  are  discussed  in  detail  in  [HBG70]. 
The  aircraft  possesses  two  sets  of  longitudinal  control  surfaces:  elevators  and 
canards.  There  are  two  available  measurements:  angle  of  attack  and  pitch  rate. 
For  our  discussion,  we  use  the  linearised  model  of  the  longitudinal  dynamics  at 
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Mach  0.0,  altitude  25,000  ft: 


x  =  Ax  +  Bu 

(5.4.1) 

-18.9  -32.1' 

-.083  0 

-2.63  0 

1.0  0  . 


[0  0  0  57.3  J 

The  plant  model  in  this  flight  condition  ia  stable  and  minimum  phase.  The 
singular  values  of  the  plant  an  shown  in  Figure  5.6.  Since  the  plant  condition 
number  is  equal  to  the  ratio  of  the  larger  to  smaller  singular  values,  we  see 
from  Figure  5.6  that  the  condition  number  is  large  at  almost  all  frequencies. 
The  physical  source  for  the  large  condition  number  is  the  ability  of  the  control 
surfaces  to  impart  relatively  large  amounts  of  rotational  energy  to  the  aircraft, 
combined  with  an  inability  to  transmit  relatively  large  amounts  of  kinetic  energy 
to  the  aircraft.  These  facts  imply  that  the  transfer  function  from  inputs  to  pitch 
rate  has  larger  gain  than  that  from  inputs  to  vertical  velocity,  or  angle  of  attack. 

The  weighting  function  for  robust  stability  (5.2.4)  fa 

fl(«)  =  i(.01e  +  l)  (5.4.2) 


y  =  Cx 


where 


A  = 


-.0226 

0 

.0123 

0 


B  = 


n  = 


-36.6 

-1.0 

-11.7 

0 


0  0 

-.414  0 

-77.8  22.4 
0  0 

0  57.3  0  0 
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while  that  for  tracking  performance  » 

'•<*>  = 


(5.4.3) 


Bode  plots  of  these  weighting  functions  appear  in  Figure  5.7,  and  a  Bode  plot  of 
the  weighted  condition  number  (5.2.19)  is  in  Figure  5.8.  As  shown  by  [Ste84], 
the  optimal  value  of  p[M]  for  this  design  is  n[M(j)u]  =  1,  Vw. 

We  shall  initially  use  weighting  obtained  by  approximating  those  given  in 
Theorem  5.3.4,  and  then  iterating  to  achieve  a  final  satisfactory  design: 


**>-[5  «£n][-»  :5][i 

0  l 

Xh-u)  1  • 
#+m  J 

(5.4.4) 

iViW  =  <-.(*)/ 

(5.4.5) 

H'jW  =  (1/S.756)r,(«)/  . 

(5.4.6) 

The  constant  factor  in  (5.4.6)  arises  from  a  crude  low  frequency  approximation 
to  the  function  /9(s)rs(s).  The  structured  singular  value  resulting  from  these 
weightings  is  plotted  in  Figure  5.9.  The  key  features  of  this  plot  are  the  two 
peaks,  one  in  the  crossover  region  and  one  at  high  frequencies  that  may  be 
reduced  by  penalising  the  control  signal  y\  more  heavily  by  adding  a  lead  filter 
to  W\: 

H',w=2-5(r?^)ri(‘)/  (s-4-7) 

while  leaving  the  other  design  weights  unchanged.  The  structured  singular  value 
of  the  second  design  (Figure  5.10)  has  a  single  large  peak  in  the  crossover  region. 
This  may  be  reduced  by  using  a  better  approximation  to  the  function  P[»)ti{»). 
Such  an  approximation  requires  that  we  introduce  a  lead  filter  into  (5.4.6): 
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Figure  5.10:  Structured  Singular  Value,  Design  2 
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(5.4.8) 


The  structured  singular  value  for  the  third  design,  using  (5.4.4),  (5.4.7),  and 
(5.4.8)  is  shown  in  Figure  5.11.  Note  that  this  plot  is  within  10%  of  the  desired 
value,  and  of  the  optimal  value  obtained  in  [Ste84].  If  desired,  one  could  reduce 
this  peak  further  through  additional  applications  of  the  techniques  developed 
in  this  chapter.  The  resulting  changes  in  the  compensator  will  only  be  minor, 
however,  and  are  likely  to  be  dominated  by  effects  not  incorporated  into  the 
linearized  design  problem. 
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Figure  5.11:  Structured  Singular  Value,  Design  3 
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Chapter  6 

Design  for  Ill-Conditioned  Plants 
Using  LQG/LTR  Method 

6.1  Introduction 

Multivariable  extensions  of  classical  loop-shaping  ideas  have  recently  received 
much  attention  [D0S8I],  [StA87],  [SLH81],  [HuM82].  The  development  of  these 
multivariable  loop-shaping  design  methodologies  was  based  upon  the  introduc¬ 
tion  of  singular  value  analysis.  Since  the  singular  value  tests  used  to  analyse 
stability  robustness  are  applicable  only  to  limited  classes  of  plant  uncertainty 
[FLC82],  [Doy82],  [DWS82],  these  design  methodologies  are  most  successful  only 
when  the  design  specifications  can  be  precisely  formulated  in  terms  of  loop  trans¬ 
fer  functions  at  one  loop  breaking  point  or  when  the  plant  has  uniform  gun  in 
all  loops. 

Owing  to  the  coupling  and  the  directionality,  feedback  properties  of  a  multi- 
variable  system  such  as  stability  margin  and  sensitivity  must  be  evaluated  at  more 
than  one  loop-breaking  point,  e.g.,  at  the  plant  input  and  plant  output.  When  the 
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plant  transfer  function,  evaluated  at  frequencies  of  interest,  is  an  ill-conditioned 
matrix,  it  can  happen  that  nominal  properties  at  one  loop-breaking  point  are  un¬ 
acceptable  even  though  properties  at  the  other  are  good.  It  can  also  happen  that 
nominal  properties  at  one  loop-breaking  point  are  extremely  sensitive  to  uncer¬ 
tainty  elsewhere  [Fre89a],  [Pre89b],  [NeM87],  [SMD88],  [Ste85j.  Hence,  when  the 
plant  is  ill-conditioned,  and  when  the  uncertainties  are  present  simultaneously  at 
different  points  in  the  feedback  loop,  multivariable  loop-shaping  techniques  such 
as  LQG/LTR  [StA87]  may  fail  to  yield  a  robust  design  [Ste84],  [Ste85]. 

Recently,  some  progress  has  been  made  in  analysing  how  to  shape  the  loop- 
transfer  functions  at  one  loop  breaking  point  to  satisfy  the  design  goals  for¬ 
mulated  in  terms  of  loop  transfer  functions  at  different  loop  breaking  points 
[Fre90]  It  is  shown  that  the  difference  between  feedback  properties  at  different 
loop  breaking  points  is  closely  related  to  the  singular  subspace  structures  of  the 
loop-transfer  functions  and  the  plant.  These  analyses  yield  conditions  that  must 
be  satisfied  by  feedback  properties  at  one  loop-breaking  point  to  prevent  poor 
feedback  properties  at  the  other  or  to  be  robust  against  uncertainty  at  the  other. 
It  turns  out  that,  to  satisfy  the  design  goals,  one  has  to  manipulate  both  the  sin¬ 
gular  values  and  the  singular  subspace  structures  of  the  loop-transfer  functions. 
However,  to  the  best  of  our  knowledge,  no  systematic  synthesis  procedures  are 
available  to  obtain  compensators  so  that  the  resultant  loops  satisfy  the  loop 
shaping  conditions  derived  in  [Fre90]. 

The  main  objective  of  this  chapter  is  to  try  to  develop  such  a  synthesis  pro¬ 
cedure.  We  shall  do  so  by  incorporating  design  insights  obtained  in  [Fre90]  into 
the  shaping  of  the  state  feedback  loop  in  the  LQG/LTR  design.  We  will  show 
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that  by  properly  selecting  weighting  functions,  the  LQG/LTR  methodology  can 
be  an  effective  loop-shaping  technique  for  ill-conditioned  plants. 

The  rest  of  the  chapter  is  organised  as  follows.  In  Section  6.2,  we  introduce 
some  notations  and  preliminaries  and  formulate  the  problem  considered  in  this 
chapter.  In  Section  6.3,  a  property  of  the  Kalman  equality  is  discussed  which 
relates  the  right  singular  vectors  of  the  return  difference  function  of  the  state 
feedback  loop  with  those  of  the  weighting  function.  In  Section  6.4,  procedures 
for  selecting  weighting  functions  are  suggested  for  the  LQG/LTR  target  state 
feedback  loop  design  which  ensure  that  the  resultant  LQG/LTR  design  has  good 
stability  robustness  against  simultaneous  uncertainties  at  the  plant  input  and 
output.  Design  examples  are  given  in  Section  6.5  to  illustrate  the  procedures. 

6.2  Preliminaries  and  Problem  Formulations 

Consider  the  linear  time-invariant  feedback  system  shown  in  Figure  6.1.  P(s) 
and  F(a)  denote  the  transfer  functions  of  the  plant  model  and  the  compensator, 
respectively.  We  assume  that  P(«)  is  an  m  x  m  rational  matrix  and  is  invertible. 
The  signals  u(s)  and  y(«)  are  the  plant  input  and  measured  output,  respec¬ 
tively.  Let  the  singular  value  decomposition  of  the  plant  model  be  denoted1 
P  =  WTZH ,  where  T  =  dtap[n,ri,...,rm]  contains  the  singular  values  with 
usual  ordering  fr[P]  =  T\  >  r*  >  •••  >  rm  =  g[P],  W  =  [u/i,u/j,. . .,u/„»]  and 
Z  =  [si.ss,...,**]  are  unitary  matrices  whose  columns  are  the  left  and  right 
singular  vectors,  respectively.  We  shall  use  a  boldface  letter  to  denote  the  col- 
1We  suppress  dependence  upon  frequency  whenever  appropriate. 


102 


iimn  space  generated  by  the  matrix  (or  vector)  denoted  by  the  corresponding 
uppercase  (lowercase)  letter.  The  condition  number  of  the  plant  is  defined  as 
<c[P]  :=  o\P)/g\P).  If  the  condition  number  is  very  large  at  some  frequency, 
we  say  the  plant  is  ill-conditioned  at  that  frequency.  The  plant  condition  num¬ 
ber  is  not  invariant  under  scaling,  or  a  change  of  units.  Hence  we  must  assume 
throughout  this  chapter  that  physically  meaningful  units  have  been  chosen.  It 
follows  that,  if  a  plant  is  ill-conditioned  with  respect  to  these  units,  then  the 
ill-conditioning  cannot  be  removed  simply  by  scaling. 

For  this  system,  the  feedback  system  properties  are  governed  by  the  following 
two  sets  of  transfer  functions: 

1.  Input  open-loop  transfer  function  Lj(a)  :=  F[a)P{a) 

Input  sensitivity  function  Sj(«)  :=  [I  +  Lj(a)]-1 

Input  complementary  sensitivity  function  2>(e)  :=  L/(a)[J  +  Li[a)]~l 

2.  Output  open-loop  transfer  function  L0  (a)  :=  P(a)F(a) 

Output  sensitivity  function  So(«)  :=  [/  +  Z/o(«)]-1 

Output  complementary  sensitivity  function  To  (a)  :=  Lo{a)\I  +  L0{a)\~l. 

When  all  design  specifications  can  be  accurately  formulated  in  terms  of  the  trans¬ 
fer  functions  at  one  loop  breaking  point,  then  multivariable  loop-shaping  design 
methods  (e.g.,  [StA87])  can  be  applied  to  manipulate  the  shape  of  the  transfer 
functions  so  that  the  design  specifications  are  satisfied  as  well  as  possible.  How¬ 
ever,  when  design  specifications  are  formulated  in  terms  of  transfer  functions  at 
different  loop  breaking  points,  then  the  multivariable  loop-shaping  methods  may 
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not  be  so  readily  applicable.  For  the  problem  we  consider  in  this  chapter,  the 
design  specifications  have  to  be  formulated  in  terms  of  transfer  functions  at  dif¬ 
ferent  loop-breaking  points.  The  following  is  the  problem  description.  It  is  the 
same  as  that  discussed  in  Chapter  5,  and  is  restated  here  for  ease  of  discussion. 

Suppose  that  the  true  plant  is  given  by 

P(>)  =  |/  +  r,(*)A,  M1P(.)[J  +  r,(»)A,«]  (6.2.1) 

where  the  unstructured  uncertainties  A<,*  =  1,2  lie  in  the  set 

D,  :=  {A  €  CTxm  :  A  is  stable  and  *[A]  <  (6.2.2) 

and  the  weighting  functions  ri(s)  and  r*(s)  describe  how  the  sise  of  each  uncer¬ 
tainty  varies  with  frequency.  Ai(s)  is  called  input  uncertainty  and  is  assumed  to 
arise  from  plant  or  actuator  modelling  errors.  A*  (a)  is  called  output  uncertainty. 
It  may  arise  either  from  modelling  error  or  from  the  use  of  fictitious  uncertainty 
to  represent  an  equivalent  performance  specification  [DWS82].  Hence,  the  prob¬ 
lem  we  are  to  consider  next  can  be  interpreted  either  as  robust  stability  against 
simultaneous  plant  input  and  output  uncertainties  or  as  robust  performance  at 
the  plant  output  against  uncertainty  at  the  plant  input.  Our  goal  is  to  find  a 
compensator  F(«),  for  given  design  specifications  ri(«)  and  rs(«),  such  that  the 
feedback  system  in  Figure  6.1  is  stable  when  the  plant  model  P(»)  is  replaced 
by  the  true  plant  P(a)  described  in  (6.1)  and  (6.2),  i.e.,  we  would  like  to  achieve 
stability  robustness  against  simultaneous  input  and  output  uncertainty.  This 
problem  is  equivalent  to  finding  a  compensator  such  that  the  structured  singular 
value  satisfies 

p[M[ju)]  <  l/Tf,Vw  (6.2.3) 
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where 


(6.2.4) 


riTt  —r\TjP~l 
—rtSoP  fj  Sq 


and  the  structured  singular  value  is  defined  at  each  frequency  by 


a 


0,  if  no  A  €  Dn  exists  with  det[/  +  M A]  =  0 
1/  min{ff[A] :  det[/  +  Af  A]  =  0  and  A  €  D«,},  otherwise. 


(6.2.5) 


However,  the  structured  singular  value  can  only  be  used  to  judge  our  compensator 
design  and  it  does  not  give  insights  into  how  we  should  construct  compensators 
such  that  (6.2.3)  is  satisfied.  To  obtain  such  insights,  bounds  on  the  structured 
singular  value  have  been  derived  (see  Chapters  2-4)  in  terms  of  the  singular 
values  and  vectors  of  the  plant  transfer  functions  together  with  the  sensitivity 
and  complementary  sensitivity  functions  used  to  perform  singular  value  analysis. 
Specifically,  we  have  for  M  defined  in  (6.2.4)  that 


max{/*i,  Ma}  <  m(M]  <  Mi  +Ma  (6.2.6) 

where 

Mi  :=  max{tf[rx7>],tf[r,S0j}  (6.2.7) 

Ma  :=  \Jo[riPSi)o[riTjP-'1}  (6.2.8) 

Hence,  to  satisfy  (6.2.3),  it  is  necessary  and  sufficient  to  keep  the  parameters  mi 
and  Ma  small. 


It  is  useful  to  decompose  an  ill-conditioned  plant  into  higher  and  lower  gain 
subsystems  at  frequencies  of  interest.  Hence  we  partition  the  plant  singular  value 
decomposition  as 

P  =  W\T\Zf  +  WxTtZf  (6.2.9) 

aFor  details  about  ■  tinctured  singular  value,  see  (Doy82j. 
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where  Tx  =  diag[TU . . . ,n),Tt  =  diag[n+u *nd  W{  and  3r,t  =  1,2,  are 
p&rtitons  of  W  and  Z  of  dimensions  compatible  with  those  of  T\  and  T%.  We 
shall  assume  that  n  >  n+l.  In  this  case,  it  can  happen  that  n%  is  very  large  even 
if  px  is  small,  which  means  that  the  system  stability  is  much  more  sensitive  to 
simultaneous  input  and  output  uncertainties  than  to  either  uncertainty  acting 
alone  [Fre89a].  Assume  that  the  gain  in  each  subsystem  given  fay  partition  (6.9) 
is  uniform,  i.e.,  fr[T<]  w  £[T<],t  =  1,2.  Then  it  is  shown  [Fre90]  that  to  prevent 
from  being  large,  the  compensator  has  to  be  chosen  so  that  the  cloaed-loop 
transfer  functions  satisfy 

9[Z?Si)9[TiZi]  C  1  (6.2.10) 

at  frequencies  for  which  the  weighting  product  |  r^j  |  is  large. 

Hence,  if  we  apply  a  multivariable  loop-shaping  design  method  at  the  plant 
input  and  would  like  to  have  stability  robustness  against  simultaneous  input 
and  output  uncertainties  or  robust  performance  at  the  plant  output,  we  must 
manipulate  the  loop  shape  to  satisfy  the  condition  (6.2.10)  as  well  as  possible  in 
addition  to  meeting  the  design  specifications  at  the  plant  input. 

6.3  Properties  of  the  LQ  Regulators 

The  LQG/LTR  procedure  consists  of  two  basic  steps.  The  first  step  is  called  the 
target  state  feedback  loop  design.  In  this  step  a  state  feedback  loop  is  designed 
to  meet  all  stability  and  performance  requirements.  The  second  step  is  called 
the  recovery  procedure  where  an  LQG/LTR  compensator  is  designed  to  asymp¬ 
totically  recover  the  state  feedback  loop.  The  quality  of  recovery  is  governed  by 
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the  nonminimum  phase  characteristics  of  the  plant  [ZhFQO].  Hence,  to  find  an 
LQG/LTR  compensator  so  that  condition  (6.2.10)  is  satisfied,  it  is  necessary  for 
the  target  state  feedback  loop  to  satisfy  condition  (6.2.10).  In  the  following  we 
shall  re-examine  some  properties  of  the  LQ  regulators  and  use  those  properties 
to  guide  us  in  the  LQ  weighting  selection  so  that  the  resulting  state  feedback 
loop  satisfies  condition  (6.2.10). 

Consider  a  state-space  realisation  of  the  plant  model  P(a) 

z(t)  =  Ax(t)  +  Bu[t)  (6.3.1) 

y(t)  =  Cx(t)  (6.3.2) 

where  x(<)  €  A”  is  the  state,  u(t)  €  R"  is  the  control  input,  y(t)  €  JF*  is 

the  measured  output.  It  is  assumed  that  (C, A)  is  observable  and  (A,B)  is 

controllable. 

Suppose  an  LQ  regulator  is  designed  for  the  system.  Hence  the  feedback  law 
is  given  by 

u(t)  =  —K,x(t)  (6.3.3) 

where  Ke  =  R~XBTP ,  and  P  is  the  unique  positive  definite  solution  of  the  Riccati 
equation 

PA  +  AtP  -  PBR~lBTP  +  BtH  =  0  (6.3.4) 

In  equation  (6.3.4)  H  is  a  design  parameter  and  is  chosen  such  that  (H,  A)  is 
observable  and  R  >  0  is  a  given  weighting  matrix. 

Denote  the  transfer  function  of  the  LQ-optimal  state  feedback  loop  by 

L(s)  4  Ka*(s)B  (6.3.5) 
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where  #(a)  =  {si  —  A)-1,  and  let  the  corresponding  aesnitivity  and  complemen¬ 
tary  sensitivity  functions  be  denoted  by 

S{»)  =  (I  +  L{s))~l  (6.3.6) 

and 

T(a)  *  £(«){/  +  £(a)j"x  (6.3.7) 

respectively.  It  is  well  known  that  L(s)  satisfies  the  following  Kalman  equality: 
[/  +  L{ju>))aR\I  +  L( /«)]  =  R  +  [fi*(jw)Blf[jr*(jw)B]  (6.3.8) 

Define  the  transfer  function  from  control  inputs  to  weighted  states  by 

0(a)  &  H{sl  -  A)~lB  (6.3.9) 

By  an  abuse  of  terminology  we  shall  call  0(a)  a  weighting  function;  our  motivar 
tion  for  this  nomenclature  is  that  we  shall  shape  properties  of  the  optimal  state 
feedback  loop  (6.3.5)  by  manipulating  0(a).  Denote  the  singular  values  of  the 
matrix  M  €  fl"*"  by  o<[M];  i  = 

Lemma  6.3.1:  Suppose  that  R  =  pi  in  (6.3.8).  Then 

(i)  Oi[I  +  L(j'w)]  =  sjl  +  {1/ p)of[Q{ju)\  (6.3.10) 

(ii)  If  Uj  is  a  right  singular  vector  of  0(j'w)  associated  with  singular  value  o<[0(ju;)], 
the  ut  is  also  a  right  singular  vector  of  J  +  L{jw)  associated  with  the  singular 
value  Oi\I  +  L{jw)\. 

Proof:  Straightforward.  ■ 
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The  first  part  of  this  lemma  is  well-known  (e.g.  [D0S8I])  and  has  been  used 
extensively  in  multivariable  loop-shaping  via  the  LQG/LTR  methodology.  The 
second  part  appears  to  be  less  well  known.  We  shall  use  this  result  to  shape  the 
right  singular  subspaces  of  the  return  difference  function  I  +  L(ju) ,  and  thus  the 
left  singular  subspaces  of  the  sensitivity  function  (6.3.6). 

Now  suppose  that  the  singular  value  decomposition  of  the  weighting  function 
(6.3.9)  may  be  separated  into  high-  and  low-gain  parts  dimensioned  compatibly 
with  those  of  the  plant  (6.2.9): 

Q  =  W.  =  [  V,,  V*  ]  [  £  ]  [  J|  ]  (6.3.11) 

where  =  dtag[o1,  •  •  •  ,oj],  =  diag[oi+u  •  •  •  ,om],  and  and  U#  are  parti¬ 

tions  of  Vq  and  Uq  whose  dimension  are  compatible  with  those  of  E,i  and  EfS. 

Proposition  6.3.2:  Assume  that  R  =  pi,  and  suppose  that  the  weighting 
function  (6.3.11)  satisfies 

efinVy/p  > 1  (6.3.12) 

and 

Uti  =  Zi ,  «*lt2  (6.3.13) 

Then  the  state  feedback  sensitivity  function  (6.3.6)  satisfies 

9[ZfS]  C  1  (6.3.14) 


Proof:  It  follows  from  (6.3.11)  and  Lemma  6.3.1  that  the  singular  value  de- 
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(6.3.15) 


composition  of  the  return  difference  matrix  can  be  written  as 


and  V  =  [  V\  Vi  ]  is  a  unitary  matrix.  Hence  the  singular  value  decomposition 
of  (6.3.6)  is 

S  =  (!  +  £)-*  =  [z,  z,  ]  [  S°-,][{£]  (6.3.16) 

The  result  follows  from  (6.3.12).  ■ 

Recall  that  the  complementary  sensitivity  function  of  the  LQ  state  feedback 
loop  must  satisfy  o\T{ju)\  <  2,  Vu>  (e.g.  [DoS81]).  It  follows  from  this  fact  and 
(6.3.14)  that  if  the  plant  is  minimum  phase  so  that  the  LTR  procedure  may  be 
successfully  applied,  then  we  can  insure  that  (6.2.10)  holds  by  requiring  that 
(6.3.12-13)  hold. 

It  is  not  possible  to  require  that  (6.3.14)  hold  at  all  frequencies,  because 
this  condition  requires  that  some  open-loop  gains  be  large.  In  particular,  it  is 
not  possible  for  loop  guns  to  be  large  at  high  frequencies.  In  such  situations, 
an  alternate  approach  is  to  require  that  &[TjZi]  <  1.  (Recall  that  fr[S(jw)]  < 
1,  Vcj.)  This  may  be  accomplished,  for  a  minimum  phase  weighting  Q(a),  by  using 
cheap  state  feedback  control.  The  following  result  is  well-known  (e.g.  [DoS81]). 

Lemma  6.3.3:  Let  R  =  pi  in  (6.3.8),  and  suppose  that  Q(s)  in  (6.3.9)  is 
minimum  phase.  Then  as  p  -*  0,  the  state  feedback  loop  transfer  function 
satisfies,  at  each  frequency, 

HM  — *  (1  (6.3.17) 
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when  W  is  a  unitary  matrix. 


Thin  result  shows  that  asymptotically  the  open-loop  state  feedback  singular 
values  satisfy  Oi[L(ju)\  «  0i[Q{jw)\/y/p  and  the  right  singular  subspaces  of  L  are 
approximately  equal  to  those  of  Q.  Suppose  that  Q  has  the  form  (6.3.11),  where 
U,i  =  Zi,i  =  1,2,  and  d[Lqi(ju)\/ ^/p  <  1.  Then  it  follows  from  the  results 
of  [FrL86],  [EreOO]  that  the  state  feedback  complementary  sensitivity  function 
satisfies  9[T(jw)Zi\  <  1.  Again,  if  the  plant  is  minimum  phase,  we  can  apply 
the  LTR  procedure  so  that  the  (6.2.10)  holds. 

Hence,  by  properly  choosing  the  weighting  Q(s)  in  the  state  feedback  loop 
design,  the  LQG/LTR  procedure  can  be  applied  to  achieve  stability  robustness 
against  simultaneous  input  and  output  uncertainties. 

6.4  Weighting  Selection  in  the  LQG/LTR  Pro¬ 
cedure 

We  can  see  from  previous  discussions  that  one  way  to  meet  requirement  (6.2.10) 
in  the  LQG/LTR  design  is  to  choose  the  LQ  weighting  function  Q(a)  so  that  it 
has  the  same  spread  in  gains  and  the  same  corresponding  right  singular  subspaces 
as  the  plant. 

It  is  clear  that  if  we  let  H  =  C,  then  Q(a)  =  C(«J  -  A)~lB  =  P(s)  and  the 
above  requirements  on  Q(a)  are  trivially  satisfied.  This  is  the  simplest  way  to 
achieve  the  requirements.  However,  it  restricts  the  loop  shapes  we  can  have  at 
the  input  loop-breaking  point.  To  have  more  freedom  in  selecting  the  loop  shape 
at  plant  input,  we  have  to  augment  proper  dynamics  to  the  plant. 
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1.  Dynamics  Augmentation  at  the  Plant  Output 

One  way  to  augment  the  plant  is  to  add  extra  dynamics  at  the  plant  output 

(see  Figure  6.2).  In  this  case  the  augmented  system  can  be  described  by 

»  «  • 

1  =A.1  +B.u,i  =  C.  *  (0.4.1) 

where 

Aa  ~  BWC  An  »  B*  ~  [  0  *  c*  =  [  0  c»  ]  (6*4*2) 

and  W(«)  :=  Cm(al  —  Aw)~lBw  is  the  dynamics  we  add  to  the  plant.  If  we 

select  H  =  [0  C*],  then  the  singular  value  and  singular  vector  properties  of  the 
resultant  LQ  feedback  loop  are  essentially  determined  by  the  augmented  plant. 
This  can  be  verified  by  writing  down  the  Equation  (6.3.8)  for  the  augmented 
system  with  the  above  weighting  and  using  Lemma  6.3.1.  We  can  see  in  this  case 
that  <?(*)  =  H(al  -  A*)~lBa  —  W(a)P(a).  Hence,  W(a)  should  be  chosen  such 
that  the  product  W(s)P(s)  preserves  the  right  singular  subspace  structure  and 
the  spread  in  gains  of  the  plant  P(a)  and  such  that  W(s)P(s)  has  the  desired 
loop  shape.  The  following  are  some  suggestions  in  how  to  choose  W[»). 

(1)  The  plant  P(s)  is  ill-conditioned  at  all  frequencies,  i.e.,  the  partition  (6.2.9) 
holds  at  all  frequencies.  In  this  case,  we  may  choose  W(s)  to  be 

W(s)  =  g(»)I 

where  g(a)  is  a  stable  minimum  phase  transfer  function  that  can  be  used  to 
manipulate  the  loop  shape  at  the  plant  input.  If  the  matrix  of  plant  left  singular 
vectors  W  can  be  well  approximated  by  a  constant  real  matrix,  8  say  W ,  then 

’When  W  wm  aufficiently  aiowly  with  frequency,  it  can  uauaUy  be  well  approximated  by  a 
real  matrix.  Technique*  of  approximating  a  complex  matrix  by  a  real  matrix  can  be  found  in 
[MaK77],  [EdK79j. 
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P(S) 


y(s) 


H  w  (s) 


Figure  6.2:  Augmenting  the  Plant  at  the  Output 


we  may  let  W  (a)  be  of  the  following  form. 


W  (a)  =  d»ag[gi(s),gt{»)t. . . ,  fc»(«)]WrT 

where  scalar  transfer  functions  ft(a),  *  =  1,2, ...  ,m,  are  selected  to  be  stable  and 
minimum  phase.  They  can  be  manipulated  to  achieve  their  desired  loop  shape. 

(2)  The  plant  P(s)  is  ill-conditioned  only  around  some  frequency  uq,  i.e.t  the 
partition  (6.2.0)  holds  only  around  frequency  wo<  A  procedure  to  obtain  the 
desired  augmenting  dynamics  is  as  follows: 

(i)  Compute  the  singular  value  decomposition  of  P{juo)  to  obtain  its  left  sin¬ 
gular  matrix  W. 

(ii)  Approximating  W  by  a  real  matrix  W. 

(iii)  Select  scalar  transfer  functions  ft(a),  »  =  l,2,...,m,  such  that  they  are 
stable  and  minimum  phase.  They  can  be  manipulated  to  achieve  the  desired 
loop  shape  at  s  =  juf o. 

(iv)  Form  the  augmenting  dynamics  W(s)  as  follows 

W{s)  =  diaglgi{8),gt{*),...,gm{s)}W 

(v)  Manipulate  free  parameters  in  ft  (a)  to  achieve  desired  loop  shape  for  W  (a)P(a). 

The  weighting  selection  procedure  will  be  repeated  in  the  LQG/LTR  design  until 
a  satisfactory  result  is  obtained. 

2.  Formal  Loop  Shaping 
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This  dynamics  augmentation  procedure  is  adopted  from  [StA87],  [Enn84] 
which  makes  it  very  easy  to  manipulate  the  loop  shape  at  one  loop  breaking 
point.  Suppose  that  the  plant  is  stable  and  the  desired  loop  shape  at  the  plant 
input  can  be  denoted  by  W(s)  =  Cw(al  —  Am)~lBw.  If  the  plant  is  augmented 
as  in  Figure  6.3,  then  the  state  space  representation  of  the  augmented  plant  is 
given  by 


A<  = 


A  0 
0  Aw 


,  B.= 


B 

Bm 


,c.=  [c  o] 


(6.4.3) 


The  LQ  optimal  feedback  gain  for  the  augmented  system  is  obtained  by  solving 
the  Riccati  equation  (6.3.4)  with  the  weighting  matrices  being  chosen  as  R  =  I 
and  H  =  [0  Cv\. 


As  we  have  shown  in  the  previous  sections,  feedback  properties  of  the  re¬ 
sultant  state  feedback  system  are  basically  determined  by  the  weighting  Q(«)  = 
H{aI—A+)~xB «  =  Cm{aI-Am)~lBw.  Note  that  by  selecting  different  Am,  Bm,Cw, 
we  now  have  complete  freedom  in  specifying  the  weighting  function  <?(«).  To  il¬ 
lustrate  how  to  select  Q(«),  we  again  consider  the  following  two  cases. 

(l)  The  plant  P(s)  is  ill-conditioned  at  all  frequencies.  In  this  case  we  may  select 


Q(«)  =  t{s)P{a) 

The  scalar  transfer  function  g(«)  can  be  manipulated  to  obtain  satisfactory  loop 
shape  at  the  plant  input.  Such  a  selection  trivially  implies  that  high  gain/low 
gain  structure  of  P  is  present  in  Q. 

(3)  The  plant  P(«)  is  ill-conditioned  only  around  some  frequency  c^.  In  this 
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case,  an  algorithm  to  obtain  the  desired  weighting  function  is  as  follows: 

(i)  Compute  the  singular  value  decomposition  of  P{ju o)  to  obtain  its  right 
singular  matrix  Z. 

m 

(ii)  Approximate  Z  by  a  real  matrix  Z. 

(iii)  Select  stable  transfer  functions  *  =  1,2,. . . ,m,  so  that  |  gi{ju)  | 
resembles  the  desired  loop  shapes  at  the  plant  input,  |  gt{jw)  |>|  0*(/u>)  I* 
if »  >  k  and  |  gi(ju0)  |  /  |  ft+i(M)  | »  /g[P(j'wo)]. 

(hr)  Form  the  weighting  function  Q(a)  as  follows 

g(«)  =  d»<v[ai(«),p»(«), . .  .*gm{»)]ZE 

Again,  the  weighting  selection  procedure  will  be  repeated  in  the  LQG/LTR  de¬ 
sign  until  a  satisfactory  result  is  obtained.  Comparing  with  the  output  dynamics 
augmentation,  the  formal  loop  shaping  method  gives  us  mare  flexibility  in  mar 
nipulating  the  loop  shapes  to  achieve  design  specifications.  However,  it  usually 
results  in  higher  order  controllers.  As  regards  to  which  dynamics  augmentation 
procedure  should  be  followed,  the  selection  may  be  based  on  whether  the  right 
or  the  left  singular  vectors  can  be  better  approximated  by  real  vectors. 

6.5  Design  Examples 

In  this  section,  we  illustrate  the  results  of  previous  sections  by  designing  a 
LQG/LTR  compensator  for  two  aircraft  control  problems. 
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Example  6.5.1  We  first  consider  the  aircraft  control  problem  studied  in  [Ste84]. 
The  plant  has  two  inputs  and  two  outputs,  and  it  is  stable  and  minimum  phase. 


It  has  strong  directionality  properties  that  vary  with  frequency.  The  problem 
formulation  is  the  same  as  the  one  in  Section  6.2.  The  state  space  representation 


of  the  plant  is  given  by 

-0.0226  -36.6000  -18.0000  -32.1000 ' 
0.0000  -1.0000  0.0830  0.0000 

A  0.0123  -11.7000  -2.6300  0.0000 

0.0000  0.0000  1.0000  0.0000 


0.0000 

0.0000 

-0.4140 

0.0000 

f  0.0000 

57.3000 

0.0000 

0.0000 

-77.8000 

22.4000 

’  °  “  [  0.0000 

0.0000 

0.0000 

57.3000 

0.0000 

0.0000  . 

The  uncertainty  is  as  described  in  (6.2.1)  and  (6.2.2)  with  weighting  functions 


given  by 


ri(«)  = 
r,(a)  = 


s  + 100 

200 
s  +  3 

2(7+003) 


Plots  of  the  plant  singular  values  and  condition  number  are  given  in  Figure  6.4- 
5.  Note  that  the  plant  is  ill-conditioned.  Plots  of  weighting  functions  r%  (a)  and 
rj(a)  are  given  in  Figure  6.6. 


Our  objective  is  to  design  an  LQG/LTR  compensator  with  the  property  that 
/i[M(jw)]  <  1/7, Vw,  for  7  as  large  as  possible.  We  shall  try  to  achieve  this 
by  manipulating  the  bounds  for  p[M],  i.e.,  the  quantities  Hi  and  p*.  From 
the  discussion  of  Section  6.2  and  6.3,  it  follows  that  as  long  as  Q  shares  the  high 
gain/low  gain  structure  of  the  plant,  we  can  simply  simply  concentrate  on  making 
Hi  small.  The  weighting  selection  can  be  judged  by  the  structured  singular  value 
plots  of  the  resulting  design. 
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magnitude  (absolute) 


magnitude  (absolute) 


.22 


For  our  case,  we  first  augment  the  plant  by  a  scalar  factor 

'  ’  <  +  0.03  <  +  3000 

at  the  plant  output.  This  will  not  change  the  gain  spread  in  each  loop  and  the 
plant  condition  number  remains  the  same. 

The  target  state  feedback  loop  is  obtained  from  (6.3.3)  and  (6.3.4)  with 
A,B,C  replaced  by  Aa,Bt,Ct  and  R  —  I  and  H  —  Cm.  The  loop-transfer 
recovery  is  obtained  using  the  Kalman  filter  with  the  parameter  g*  =  1000.  The 
structured  singular  value  of  the  resulting  design  is  given  in  Figure  6.7.  One  can 
see  that  for  this  design,  n[M(ju))\  <  1.28,  Vw,  i.e.,  the  system  will  remain  stable 
for  simultaneous  input  and  output  uncertainties  whose  magnitudes  are  leas  than 
«  0.8.  Hence,  a  decent  result  can  be  obtained  throurv  relatively  simple 
LQG/LTR  procedure. 

Example  6.6.2  We  now  consider  the  drone  lateral  attitude  control  problem  stud¬ 
ied  by  Ridgely  and  Banda  [RiB86].  We  shall  use  this  example  to  illustrate  that 
when  the  plant  is  ill-conditioned,  the  singular  value  balancing  method  suggested 
in  [Ath86]  may  result  in  a  nonrobust  design.  The  state  space  representation  of 
the  plant  is  given  by 

-0.08527  -0.0001423  -0.9994  0.04142  0  0.1862 

-46.86  -2.757  0.3896  0  -124.3  128.6 

-0.4248  -0.06224  -0.06714  0  -8.792  -20.46 

0  1  0.0523  0  0  0 

0  0  0  0  -20  0 

0  0  0  0  0  -20 . 

[000020  0lT  [l  0  0  0  0  0' 

[0000  0  20J  ,u‘[o  0  0  10  0 
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magnitude  (absolute) 


The  plant  is  minimum  phase,  but  unstable.  The  design  specifications  can  be 
represented  by  the  following  weighting  functions  on  the  input  and  output  uncer¬ 
tainties 


n(«)  = 


0.5s  + 1 
10 


„(,) .  «fe±a. 


Plots  of  the  plant  singular  values  and  condition  number  are  given  in  Figures 
6.8-0.  Note  that  the  plant  is  ill-conditioned.  Plots  of  weighting  functions  rx(s) 


and  rs(s)  are  given  in  Figure  6.10. 


First,  we  examine  two  designs  given  in  [RiB86].  Design  1  (see  Table  13.6  in 
[RiB86])  uses  the  singular  value  balancing  method  where  the  plant  is  augmented 
by  adding  an  integrator  in  each  loop  and  the  weighting  function  in  the  LQ  design 
was  chosen  in  such  a  way  that  the  singular  values  of  the  state  feedback  loop  were 
squeesed  together  at  low  frequencies.  The  authors  of  [RiB86]  noticed  that  even 
though  this  design  satisfied  design  specifications  at  the  plant  input,  it  failed 
to  meet  design  specifications  at  the  plant  output.  Design  2  (see  Table  13.7  in 
[RiB86])  also  augments  the  plant  by  adding  an  integrator  in  each  loop.  In  this 
case  the  weighting  functions  in  the  LQ  design  is  simply  chosen  as  the  plant,  i.e., 
Q(s)  =  C(al  —  A)~lB  in  our  notation.  The  authors  of  [RiB86]  noticed  that,  for 
this  design,  the  design  specifications  were  not  only  satisfied  at  the  plant  input 
but  also  at  the  plant  output.  Using  our  results  in  previous  sections,  we  are  able 
to  explain  the  above  observations.  Clearly,  the  weighting  selection  in  Design  2 
happens  to  satisfy  the  directionality  requirements,  hence  the  resulting  design 
meets  condition  (6.2.10),  which  implies  a  robust  design  against  both  input  and 
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Figure  6.8:  Plant  Singular  Values,  Example  6.5.2 
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magnitude  (absolute) 


magnitude  (absolute) 


output  uncertainties.  On  the  other  hand,  since  the  weighting  selection  of  Design  1 
does  not  meet  these  directionality  conditions  and  the  plant  is  ill-conditioned,  it 
is  not  surprising  that  Design  1  results  in  a  nonrobust  design.  In  Figures  6.11-12, 
we  plot  the  bounds  for  the  structured  singular  values  of  the  two  designs.  One 
can  see  that  Design  1  indeed  has  much  larger  structured  singular  values  than 
Design  2. 


Next,  we  illustrate  that,  by  adding  a  little  more  dynamics  to  the  plant  output, 
we  can  further  improve  Design  2.  We  first  augment  the  plant  by  a  scalar  factor 

,  _  10(s  +  2)(s  +  3) 

s(«  +  200) 

at  the  plant  output.  The  target  state  feedback  loop  is  obtained  from  (6.3.3) 
and  (6.3.4)  with  A,B,C  replaced  by  AmiBmyCm  and  R  —  I  and  H  =  Cm.  This 
weighting  selection  satisfies  the  directionality  criteria.  The  loop-transfer  recovery 
is  obtained  using  the  Kalman  filter  with  the  parameter  q  =  1000.  The  bounds  for 
structured  singular  values  of  the  resulting  design  is  given  in  Figure  6.13.  One  can 
see  that  for  this  design,  the  upper  bound  for  fx[M{ju)\  is  less  than  1.66,  Vw,  while 
the  lower  bound  for  n[M(jw)\  of  Design  2  is  greater  than  1.7  at  some  frequencies. 


6.6  Concluding  Remarks 

In  this  chapter,  we  have  considered  the  problem  of  obtaining  stability  robustness 
in  the  presence  of  simultaneous  input  and  output  uncertainties.  We  demonstrate 
how  the  LQG/LTR  design  methodology  can  be  effectively  applied  to  this  prob¬ 
lem.  The  key  in  the  application  is  the  weighting  selection  in  the  LQ  optimal  state 
feedback  loop  design.  It  is  shown  that  the  right  singular  subspaces  of  the  return 
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Figure  6.11:  Bounds  for  SSV,  Design  1,  Example  6.5.2 
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Figure  6.12:  Bounds  for  SSV,  Design  2,  Example  6.5.2 
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Figure  6.13:  Bounds  for  SSV,  Design  3,  Example  6.5.2 
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difference  matrix  of  the  target  LQ  optimal  state  feedback  loop  are  uniquely  de¬ 
termined  by  that  of  the  weighting  functions  of  the  LQ  design.  This  provides  an 
effective  means  to  obtain  a  compensator  which  achieves  the  loop-shaping  con¬ 
ditions  derived  in  (FTe90].  To  achieve  stability  robustness  against  simultaneous 
input  and  output  uncertainties,  the  singular  subspace  structure  of  the  plant  must 
be  taken  into  account  in  the  weighting  selection.  Hence,  both  singular  values  and 
singular  subspace  structures  of  the  loop-transfer  functions  should  be  manipulated 
to  achieve  a  satisfactory  design. 
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Chapter  7 

Conclusions  and 
Recommendations  for  Future 
Research 


A  number  of  problems  remain  before  a  satisfactory  robust  multivariable  design 
methodology  is  available.  Some  of  these  problems  are  now  described. 

It  appears  that  ill-conditioned  plants  may  be  inherently  difficult  to  robustly 
control  for  some  classes  of  uncertainty.  However,  the  sise  of  the  plant  condition 
number  depends  upon  the  units  chosen  with  which  to  measure  the  plant  input 
and  output  signals.  This  fact  is  somewhat  disturbing,  as  inherent  difficulty  of 
a  design  problem  should  be  independent  of  units.  A  little  thought,  however, 
suggests  that  changing  units  to  reduce  the  sise  of  the  plant  condition  number 
may  increase  the  sise  of  the  plant  uncertainty.  Hence  the  difficulty  of  the  design 
problem  would  be  preserved  under  scaling.  Further  research  is  necessary  to 
resolve  this  matter. 

The  results  described  above  focus  on  analysing  system  properties  at  a  fixed 
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frequency.  Yet  it  is  crucial  to  understand  design  tradeoffs  that  must  be  performed 
between  system  properties  in  different  frequency  ranges.  For  example,  it  appears 
that  the  inherent  difficulty  associated  with  ill-conditioned  plants  takes  the  form 
of  a  frequency  dependent  design  tradeoff.  A  tentative  conjecture  is  that  the 
loop  singular  values  must  be  separated  by  an  amount  proportional  to  the  plant 
condition  number,  and  that  the  requirement  for  this  separation  aggravates  the 
tradeoff  imposed  by  the  Bode  gain-phase  relation.  Verifying  this  conjecture  seems 
to  require  use  of  the  multivariable  gain-phase  relations  developed  in  [FrL:88]. 

Design  limitations  that  are  due  to  unstable  poles  and  nonminimum  phase 
seros  are  fundamental  in  classical  control  design.  These  tradeoffs  are  governed 
by  integral  relations.  An  important  avenue  of  research  is  to  extend  these  results 
to  multivariable  systems  using  the  mathematical  tools  developed  in  [JVL88]. 

The  preceding  analysis  problems,  and  others,  appear  to  lie  at  the  heart  of 
understanding  multivariable  design  tradeoffs.  Once  such  an  understanding  is 
available,  it  will  be  possible  to  incorporate  it  into  design  methodologies.  This 
should  proceed  along  two  lines.  First,  for  systems  with  two  or  three  inputs 
and  outputs  it  should  be  possible  to  develop  a  direct  multivariable  loop-shaping 
methodology.  Secondly,  as  was  done  in  Chapters  5  and  6  of  this  report,  methods 
must  be  found  to  manipulate  feedback  properties  expeditiously  through  weighting 
function  selection  for  formal  synthesis  procedures. 
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Appendix  A.  Proof  of  Proposition  3.5.4 


We  shall  first  need  several  preliminary  lemmas  in  order  to  prove  Proposition  5.4.  For 

*11  *u 
*11  *» 

Lemma  A.1  Let  M  be  given  in  (9.4-9)  with  n  =  2,  k  =  1.  Then, 


convenience,  we  shall  also  denote  Tj  = 


rank[M]  <  2 


(A.1) 


*1 


Proof:  It  can  be  readily  verified  that 


M  =  diag[ru  rj,  ra/j]  • 


TjP~  i 
-So 


[P  -  /,]. 


Lemma  A.2  Let  M  be  given  in  (9.4.9)  with  n  —  2,  k  —  1,  an 4  let  V  =  diag[ «i,  Uj,  u8], 
V  =  diag[vi ,  vj,  ife],  where  ult  uj,  vXj  vj  G  C,  uSl  vs  G  C*xl,  ||u,||  =  ||v<||  =  1.  Then, 

rank\VaMU ]  <  2  (A.2) 

Proof:  Noting  that  ranfc[{7]  =  ranJb[V]  =  S,  and  ranfc[Af]  <  2,  the  lemma  follows 
immediately  from  the  well-known  Sylvester  inequality. 


Lemma  A.S  Let  the  interaction  parametere  be  given  in  (9.4.11-19),  and  (9.4-17-18). 
Sappoee  that  n  =  2,  4=1.  Then  theee  parametere  are  related  via 

(i) 

p»(l,2,S)  <  p*(l,S,2)  +  m*(1,2)  .ff[r,5o]  +  p*(2,5).  |  n*u  | 

+  p*(l,5)*  |  rjtjj  |  +  |  ri*n  |  •  |  rjtu  |  »[ra5o]  (A.S) 


4 


% 
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(«) 

/(1.S.2)  <  /i8(l,2,3)  +  MJ(l,2).»[r8S0]  +  M,(2l3).|fit11| 

+  I  I  +  |  ritn  |  •  |  fj*»  |  ff[r»5o]  (A.4) 

Proof:  By  symmetry,  it  suffices  to  show  Part  (i).  Toward  this  end,  choose  t»i  =  uj  = 
t>i  =  vi  =  1,  «a  =  (rjefP_1To)F/||rj«fP“1To||l  and  t*  =  (raSoPei)/IKsSoPei||  “ 
Lemma  A.2.  Consider  the  Laplace  expansion  of  det\yaMU]  =  0.  The  result  follows  fay 
noting  that 

I  ritu  ||  rahi  ||  vffrsSou$  |  <  /** (1, 2)  •  *[rjS0] 

|  ri*n  ||  vfrtSoPet  |  »[rje?P~lTo]  <  /ia(2,3)-  I  Mil  I 
I  ritn  ||  ri«f  P_1T0us  |  ff[rs5oPe|]  <  /i*(I,3)-  |  Mm  | 

I  Mil  ||  vf  r8S0ea  ||  nef P'^otia  |  <  /i*(l,3,2) 

I  Mil  ||  Mm  II  wfrjSouj  |  <  |  Mil  H  Mm  I  **[»lSb] 

and  applying  the  triangle  inequality. 


f 


Proof  of  Proposition  3.5.4 

Proof:  By  symmetry,  we  assume  /i(l,  2, 3)  =  /is  with  no  loss  of  generality.  Then,  it 
suffices  to  prove 

#*(1,2.3)  <  1.62  max  {m,  /ij>  (A.5) 


Noticing  that 

/i*(l,2,3)  •  /is(l,3,2)  =  #**(1,2)  •  m*(1,8)  •  #**(2.3)  (A.6) 
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Solving  this  inequality  yields 


This  completes  the  proof. 


*4<(\/8  +  2)fs 
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Appendix  B.  Proof  of  Proposition  3.5.6 


Proof:  First  note  that  the  lower  bound  in  (3.5.15)  is  a  restatement  of  (3.4.14).  In  order 

to  show  the  upper  bound,  consider  instead  that  M  =  Hag  [Mix,  Ala]  +  Mi,  where 

0  0  Mu 

Mi  =00  Mu 
Mm  Mm  0 

It  follows  from  Chapter  2  that  /i[M]  <  fix  +  m[Ml].  From  [?],  it  follows  that  fi[Mx]  = 
maxx.y  p[XaM\Y]  for  X  =  dtog[xi,  xa,  xj]  and  K  =  diag[ yj,  yj,  y»],  where  Zj  and 
p,-  are  unit  vectors  with  appropriate  dimensions.  It  can  then  be  readily  verified  that 
M[Mx]  <  >/mS(1>8)  +  M*(2.3).  The  proof  is  now  completed  by  noticing  that  (3.5.16*18) 
follow  from  straightforward  calculations. 
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Appendix  C.  Proof  of  Proposition  4.4.2 


Before  proving  Proposition  4.4.2,  we  first  consider  s  partition 


where  M(/)  satisfies 

mo = 


m=E«(0,  (c.i) 

1=1 


My  Fori=lt- +  jsi  +  f-1; 

My  Fori=*N  +  2-lt-",N,  j  =  i  +  l-l-N-,  (C.2) 

0  Otherwise 


Note  that  M(i)  has  all  blocks  being  sero,  except  Mu,  Mtli+i,  •  •  • ,  Mff+i-ijf,  Mnjl-i  •  •  •, 
Mff+i-t.i-  Specifically,  there  is  only  one  block  possibly  being  nansero  along  each  block- 
row  and  block-column.  Td  illustrate,  consider  N  =  4. 

Example  C.I:  Let  N  =  4.  Then, 


Af(l)  =  *ag[Mn,  Mn,  Mn,  Mu],  M(2)  = 


M(3)  = 


0  0  Mu  0 

0  0  0  Mu 

Mu  0  0  0 

0  Mo  0  0 


M(4)  = 


0 

Afu 

0 

0 

0 

0 

Mn 

0 

0 

0 

0 

Mu 

Mil 

0 

0 

0 

a 

0 

0 

0 

MU 

M„ 

0 

0 

0 

0 

Mn 

0 

0 

0 

0 

Mn 

0 

Proof  of  Proportion  44.B:  With  no  loss  of  generality,  we  assume  ns  =  (ffuffn  •  •  •  owl)1/*. 
Note  first  that 


=  inf  max 
*i, —4h 


U'ui  s'1* 

The  solution  to  this  minimisation  problem  is  obtained  by  setting 

«?i  =  1 

*  =  Vi  >  2. 
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This  leads  to 


i 

r 
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